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ABSTRACT 


The  multi-period  open  pit  mine  production  scheduling  problem 
is  formulated  as  a  large  scale  linear  programming  problem 
using  the  block  concept.  A  solution  procedure  is  developed 
through  decomposition  and  partitioning  of  the  subproblem  into 
elementary  profit  routing  problems  for  which  an  algorithm  is 
presented.  Many  of  the  traditional  mine  planning  concepts 
are  discussed  and  suggestions  for  improvement  through  use  of 
the  techniques  developed  in  this  thesis  are  given. 

In  the  development  of  the  solution  procedure,  those 
constraints  which  govern  the  mining  system  are  considered  as 
the  master  problem.  The  constraints  which  dictate  the 
sequence  of  extraction  are  used  as  the  subproblem.  The 
properties  of  the  single  period  subproblem  and  its  dual  are 
discussed,  and  the  dual  problem  is  shown  to  be  equivalent  to 
a  bipartite  maximum  flow  problem  for  which  an  algorithm  is 
given.  The  multi-period  subproblem  algorithm  is  developed  by 
partitioning  by  stages  and  using  the  properties  of  the  single 
period  subproblen. 

This  treatment  allows  optimization  of  the  complete  mining- 
concentrating-ref  inlng  system  over  the  entire  planning 
horizon  and  permits  the  system  to  dictate  how  and  when  to 
process  a  block  of  material. 
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CHAPTER  1 


INTRODUCTION 

1.1  MINE  PLANNING 

The  mining  of  mineral  deposits  in  such  a  manner  that  at  depletion  the 
maximum  possible  profit  is  realized  has  been  an  unsolved  problem  since  man's 
discovery  of  the  usable  elements  burled  beneath  the  earth's  surface.  In 
the  days  when  high  grade  reserves  were  adequate  to  supply  our  needs  the 
attention  given  this  problem  was  negligible.  The  philosophy  at  that  time 
was  to  extract  the  material  in  an  orderly  fashion,  keeping  in  the  high  grade 
until  depletion.  Right  or  wrong,  profits  were  high,  so  no  question  of 
optimum  profitability  confronted  the  operators.  Since  World  War  II  and  the 
depletion  of  the  most  accessible  of  the  world's  high  grade  reserves,  the 
mining  Industry  has  been  forced  into  working  with  lower  grade  material. 

The  sequence  of  extraction  has  now  become  more  Important;  and,  in  many  cases, 
has  become  a  problem  whose  solution  is  vital  to  the  existence  of  a 
profitable  operation. 

The  planning  of  an  extraction  sequence  over  a  particular  time  horizon, 
typically  the  life  of  the  mineral  deposit,  is  commonly  referred  to  as 
mine  planning.  Mine  planning  is  usually  divided  into  three  categories; 
long  range,  short  range ,  and  operational  [8]. 

A  long  range  plan  defines  the  ultimate  economic  limit,  or  optimum 
pit  limit,  i.e.,  defines  the  size  and  shape  of  an  open  pit  at  the  end  of 
its  life  [8],  [10],  [16],  [20].  The  long  range  plan  serves  as  an  aid  in 
the  evaluation  of  the  economic  potential  of  a  mineral  deposit  and  delineates 
the  economic  ore  body.  This  analysis  is  essential  in  the  planning  for 
surface  facilities  such  as  treatment  plants,  waste  dump,  tailing  ponds  and 
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other  elements  complementary  to  the  mining  operation.  In  some  instances 
the  long  range  plan  also  serves  as  a  guide  for  short  range  plans. 

Short  range  plans  are  a  sequence  of  depiction  schedules  leading  frofe 
the  initial  condition  of  the  deposit  to  the  ultimate  pit  limit  [8],  [1]. 

These  plans  arc  developed  subject  to  physical,  geological,  operating,  legal, 
and  other  policy  constraints.  Each  plan  usually  varies  in  duration  from 
one  to  ten  years  and  provides  information  necessary  for  forecasting  future 
production  and  capital  expenditures. 

Operational  or  actual  production  planning  is  concerned  with  the  present 
operating  state  within  the  confines  of  the  most  recent  short  range  plan. 

It  is  the  operator's  guide  for  orderly  mining  to  gain  the  objectives  of  the 
present  as  well  as  the  short  range  plans  under  the  constraints  of  present 
conditions  and  policies.  The  planning  period  is  usually  a  year  with  stages 
of  months,  weeks,  or  days. 

The  division  of  mine  planning  into  the  categories  described  above 
closely  follows  the  stages  of  evaluation,  data  refinement  and  development 
of  an  ore  deposit.  For  example,  in  the  initial  stages  of  an  evaluation, 
only  sufficient  data  is  obtained  from  drilling,  geological  studies  and 
other  sources,  to  determine  the  economic  feasibility  of  the  deposit  and 
hence  whether  or  not  it  seems  profitable  to  continue  to  invest  time  and 
money  in  a  particular  deposit.  If  such  an  Investment  is  deemed  profitable, 
more  data  is  obtained  through  additional  expenditures.  With  this  additional 
Information  the  long  range  plans  are  refined,  and  the  evaluation  extends 
to  the  short  range  planning  stage. 

Operational  planning  is  the  final  stage  of  planning  and  is  based  on 
the  best  possible  information  available;  that  is,  information  from  exploration, 
development  and  production  drilling,  extensive  geological  studies,  production 
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records,  economic  studies  and  forecasts  and  other  sources. 

Since  the  characteristics  of  the  material  to  be  m'.ned  arc  based 
primarily  on  a  sampling  program  and  also  since  the  mining  process  typically 
extends  over  a  considerable  number  of  years,  there  exists  a  great  deal  of 
uncertainty  in  the  physical,  economic  and  technological  factors  upon  which 
planning  is  based.  Also,  management's  objectives  may  vary  over  time  due 
to  changes  in  economics  and  technological  conditions.  Tnere.  is  considerable 
interaction  between  the  categories  of  mine  planning  each  of  which  may  have 
a  different  objective.  Clearly  as  we  progress  from  long  range  to  operational 
planning,  the  degree  of  uncertainty  decreases. 

For  the  reasons  of  uncertainty,  changing  objectives,  and  the  increasing 
availability  of  more  and  better  information  as  the  operation  progresses 
discussed  in  the  previous  paragraph,  the  necessity  of  the  different  categories 
of  mine  planning  and  their  continual  revision  is  quite  evident. 

As  additional  data  becomes  available  and  objectives  change  each  category 
of  planning  is  updated  so  as  to  reflect  the  refined  information  or  changing 
conditions  of  the  times. 

The  primary  problem  to  be  considered  in  this  investigation  is  open 
pit  mine  planning  in  the  short  range  and  operational  stages.  This  is  what 
we  will  call  open  pit  mine  production  scheduling.  As  will  become  clear  to 
the  reader,  as  he  progresses  through  the  development  of  the  approach  to  this 
problem  given  in  this  study,  the  effect  of  uncertainty  will  be  lessened; 
changes  in  economic  and  technological  conditions  and  the  integration  of  the 
planning  stages  may  be  accomplished  with  relative  ease  with  our  method. 


1.2  TRADTTJON’Al  APPROACH  TO  MINK  PLANNING 


The  traditional  method  of  mine  planning  is  the  trial-and-error , 
hand-calculated,  cross-section  approach.  This  method  consists  of  tracing 
a  trial  pit  limit  on  vertical  and/or  horizontal  sections,  taking  into 
consideration  wall  slopes,  geology  and  plan  objective.  The  material  within 
the  trial  pit  is  divided  into  ore  and  waste  on  the  basis  of  geological 
interpretations  and  an  economic  or  grade  cut-off.^  Volumes  of  material 
are  determined  by  planimetering  the  areas  on  the  sections  and  multiplying 
by  a  proper  factor  based  on  the  distance  between  sections  and  geological 
interpretations. 

Grade  analyses  are  assigned  to  the  various  types  of  material  based  on 
drill-hole  data  by  various  method  of  assigning  drill-hole  sample  Influence 
to  a  certain  volume  of  material.  The  trial  pits  are  expanded  or  contracted 
to  meet  desired  requirements  by  considering  the  profitability  of  small 
increments  surrounding  the  trial  pit. 

The  profitability  of  the  increments  and  trial  pit  is  determined  by  a 
ratio  of  cubic  yards  of  waste  to  tons  of  ore.  This  ratio  is  usually 
referred  to  as  the  break-even  stripping  ratio  [20],  or  bottom  stripping 
ratio  [10].  This  stripping  ratio  indicates  the  point  where  it  is  uneconomical 
to  remove  ore  considering  the  amount  of  waste  this  necessitates  in  removing. 
The  basis  for  this  ratio  is  usually  purely  economic  or  empirical. 

Obviously,  there  are  many  faults  in  this  method,  and  also  in  the  basic 
foundations  on  which  the  many  judgements  involved  are  made,  in  particular, 
cut-off  grade  and  stripping  ratio.  There  shortcomings  are  discussed  in 
[1],  [10],  [19],  and  [20].  The  advent  of  computers  Improved  mine  planning 
methods  considerably  as  discussed  in  the  following  section. 

+Cut-off  grade:  the  criterion,  with  an  economic  motive,  normally  evploytfd 
in  mining  to  discriminate  between  ore  and  waste  in  a  mineral  deposit  [14], 

[22] .  (Section  1.4) . 
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1.3  GENERAL  DESCRIPTION  OF  PROBLEM 
1.3.1  Problem  Area 

Many  of  the  shortcomings  of  the  traditional  long  and  short  range  mine 
planning  methods  have  been  overcome  in  the  past  few  years.  Much  of  the 
Improvement  can  be  attributed  solely  to  the  speed  of  data  manipulation  and 
calculation  with  computers,  since  many  of  the  methods,  as  reported  in  the 
literature  [1],  [20],  are  patterned  after  the  conventional  methods. 

Elimination  of  the  emphasis  on  the  stripping  ratio  is  one  of  the  improvements 
made.  The  "block  concept,"  which  will  be  discussed  in  more  detail  later, 
la  a  key  to  the  Improvements  in  the  planning  methods  [1].  Very  little  has 
been  accomplished  in  the  area  of  operational  planning. 

The  operational  mining  system  may  be  considered  as  a  combination  of 
three  subsystems;  mining,  concentrating,  and  refining.  The  mine  or  group 
of  mines  is  the  area  in  which  the  material  is  extracted  from  its  position 
in  the  earth  and  this  subsystem  also  usually  Includes  the  transportation 
to  the  concentrator  or  waste  dump.  The  concentrator  is  a  plant  or  group 
of  plants  where  the  crude  ore  is  upgraded  by  various  processes,  depending 
on  the  type  of  ore,  to  a  concentrate  which  is  amenable  for  refining.  The 
refinery  contains  the  finishing  process  which  produces  the  product  ready  for 
the  manufacturing  market.  Not  all  operations  include  all  three  subsystems 
and  some  may  even  Include  more,  such  as  the  marketing  phase.  The  model 
discussed  in  this  work  is  general  enough  so  that  it  conforms  to  any 
particular  operation. 

The  goal  of  management  is  usually  to  maximize  some  form  of  profit 
subject  to  the  particular  constraints  of  their  operation.  The  form  of  profits 
is  usually  total  profits,  present  value,  or  immediate  profit. 


The  purpose  of  our  research  Is  to  develop  an  approach  to  production 
scheduling  which  Improves  on  the  trial-ond-error  concept  still  prevalent 
in  present  techniques,  and  does  away  with  the  present  concept  of  cut-off 
grade  or  economic  cut-off  based  on  purely  marginal  analysis.  This  cut-off 
concept  will  be  discussed  in  Section  1.4  of  this  chapter. 

The  objective  of  production  scheduling  is  .o  determine  a  feasible 
extraction  schedule  which  maximizes  profits  over  the  planning  period. 

(The  term  profits  here  may  be  actually  revenue  minus  costs,  or  present  vslue; 
differentiation  of  these  values  is  not  necessary  for  purposes  of  this 
investigation.)  A  feasible  schedule  is  one  which  satisfies  a  number  of 
constraints  on  factors  such  as:  orderly  extraction,  mining  equipment 
capacity,  milling  capacity,  refining  or  market  capacity,  grades  of  mill 
feed  and  concentrates,  labor,  and  other  physical,  operating,  legal,  and 
policy  limitations. 

To  develop  a  model  to  investigate  this  problem,  we  need  to  expand 
on  a  key  factor,  the  block  concept. 


1.3.2  The  Block  Concept 

First  of  all,  consider  the  actual  methodology  of  open  pit  mining.  Due 
to  the  type  of  equipment  available  and  its  capabilities,  open  pit  mining 
usually  proceeds  along  a  series  of  benches  as  shown  in  Figure  1.1.  The 
mining  sequence  usually  consists  of  blasting  a  certain  volume  of  material  In 
a  bench,  loading  this  material  into  a  haulage  unit,  and  transporting  It  to  a 
waste  dump  or  concentrator.  The  blasting  patterns  usually  are  rectangular 
and  material  from  the  entire  height  of  the  bank  is  loaded  Into  a  unit.  Thus 
a  very  accurate  description  of  a  mining  operation  can  be  given  in  terms  of 


FIGURE  1.1:  OPEN  PIT  MINING  OPERATION 
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mine-able  units  or  "blocks" — remove  a  block  of  material  from  the  deposit 
and  transport  it  to  its  destination.  One  of  the  earliest  reports  on  the 
"block  concept"  is  given  in  [1],  Since  the  actual  mining  is  accomplished 
by  blocks,  it  follows  that  the  concentrating  and  refining  can  be  similarly 
described  as  a  treatment  of  blocks.  This  idea  is  illustrated  in  Figure  1.2. 
Note  the.  reduction  in  block  size  after  concentrator,  illustrating  the 
removal  of  some  waste  material  from  the  block. 

Since  the  process  can  be  so  accurately  described  in  terms  of  blocks, 
it  seems  clear  that  the  block  concept  should  provide  a  very  good  estimate 
of  the  input  and  output  to  the  system.  The  "block"  then  becomes  the  logical 
unit  to  which  production  scheduling  can  be  applied. 

The  block  size  is  influenced  by  equipment  capabilities,  geological 

structure,  allowable  wall  slopes,  accuracy  of  sample  data,  manner  of  mining, 
desired  use  of  block,  and  capability  of  manipulating  a  huge  number  of  blocks 
[1],  [11].  The  smaller  the  block  size,  the  more  flexibility  in  planning 
is  available  and  also  more  refinement  possible,  but  the  number  of  blocks 
Increases.  Actually,  very  small  blocks  are  not  very  realistic,  since  the 
economic  and  assay  data  upon  which  scheduling  is  dependent  are  based  on 
samples  from  drill-holes  which  may  range  in  separation  from  40  to  1000  feet. 
The  height  of  a  block  should  usually  be  taken  as  the  bench  height  since 
this  is  the  way  it  is  mined.  Examples  of  block  dimensions  that  have  been  used 
in  previous  studies  are:  100'  x  100'  x  40',  40'  x  40'  x  40'  and 
100'  x  50'  x  50  [1],  [10],  [20]. 

1.3.3  Preliminary  Block  Development 

The  manner  in  which  a  deposit  is  divided  into  blocks  is  a  function  of 
criteria  similar  to  the  block  size.  One  could  accomplish  this  by  hand  and 


FIGURE  1 . 2 
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thus  consider  a  great  deal  of  geological  detail,  but  the  more  realistic 
way  Is  to  do  the  job  with  a  computer.  The  reduction  in  labor  usually  far 
outweighs  any  refinement  gained  in  the  manual  technique. 

Once  the  deposit  has  been  divided  into  blocks,  the  geologieal,  physical, 
and  economic  characteristics  of  each  block  must  be  obtained  before  production 
scheduling  can  be  undertaken.  That  is,  a  mineralised  block  Inventory 
providing  the  following  data  must  be  developed: 

1.  Block  Identification 

2.  Volumes  (crude,  concentrate,  etc.) 

3.  Analysis 

4.  Material  Classification 

5.  Mining  Equipment  Hours 

6.  Treatment  Plant  Hours 

7.  Economics 

The  development  of  this  data  is  based  on  drill-hole  samples,  metallurgical 
test  data,  and  economic  and  technological  features  of  prdsent  or  proposed 
mining  systems.  It  is  not  the  purpose  of  this  work  to  go  into  the  development 
of  this  data  since  the  methods  are  peculiar  to  each  operation  and  examples 
of  such  methods  are  given  in  [1],  (8],  (19),  and  (20). 

Given  this  block  inventory  which  provides  the  necessary  economic, 
technological,  and  geological  data,  we  are  ready  for  evaluation  studies  and 
production  scheduling. 
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1.4  Till-.  CUT-OFF  CON’CKPT 

The  cut-off  concept  is  nn  economical ly-b.tsed  criterion  which  is 
normally  used  in  mining  to  discriminate  between  waste  and  ore  in  a  mineral 
deposit  [8],  [14],  [19].  As  presently  employed,  it  is  usually  n  static 
cut-off  as  compared  to  a  dynamic  economic  cut-off  which  will  be  proposed 
in  this  discussion.  The  essential  difference  is  that  the  dynamic  cut-off 
is  a  function  of  the  state  of  the  mining  system  at  the  time  a  decision  is  to 
be  made  as  well  as  the  future  effects  of  this  decision.  This  difference 
will  be  made  clear  in  what  follows.  The  profits  of  a  mining  operation  can 
be  appreciably  affected  by  the  choice  of  a  cut-off  and  hence  this  is  an 
important  consideration. 

The  traditional  cut-off  grade  is  a  product  of  static  or  semi-dynamic 
marginal  analysis.  An  extremely  simple  example  of  this  is  shown  graphically 
in  Figure  1.3.  In  this  figure,  A  could  be  any  parameter  with  which  one 
wants  to  judge  profits  such  as  cut-off  grade,  volume,  depth,  etc.!  If  X 
represents  cut-off  grade,  then  would  be  the  optimum  cut-off  grade. 

Note  Aq  is  the  point  where  marginal  revenue  equals  marginal  cost. 

This  analysis  is  usually  made  without  regard  to  the  state  of  the 
mining  system.  Thus  the  cut-off  grade  is  essentially  determined  independent 
of  the  mining  sequence,  capabilities  of  the  mining  system,  and  other 
operational  constraints.  K.  F.  Lane  recognized  the  fallacy  of  this  approach 
and  proposed  in  [14]  that  the  capacities  of  the  various  subsystems  should  be 
considered  as  well  as  economics  when  determining  a  cut-off  grade.  He 
assumed  the  mining  sequence  to  be  known,  even  though  the  sequence  may  be 
influenced  by  the  cut-off  choice,  since  the  tools  were  not  available  for  him 
to  discuss  the  nore  complex  problem. 

The  economic  cut-off  which  may  be  atrlbuted  to  the  block  concept  has 
replaced  the  grade  cut-off  in  some  operations  [1],  [8],  [10].  It  is  an 
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improvement  over  the  traditional  out-off  grade  in  that  it  gives  weight 
to  factors  such  as  location,  material  character 1st Jet;  which  affect  operations, 
and  other  factors  which  are  not  easily  taken  care  of  with  the  grade 
approach.  The  main  objection  to  its  use  is  that  it  is  normally  a 
prejudgement.  That  is,  it  is  not  based  on  the  influencing  factors  of  the 
mining  system.  Usually  a  volume  of  material  (block)  is  classified  as 
ore  or  waste  on  the  economics  or  characteristics  of  the  single  block  itself 
without  considering  the  system  or  the  interaction  of  other  blocks.  At  times 
this  practice  may  be  justified,  but  in  general  it  is  in  error. 

As  a  simple  example  of  where  the  static  economic  cut-off  breaks  down, 
consider  the  situation  illustrated  in  Figure  1.4.  The  value  of  the  blocks 
as  both  an  ore  and  a  waste  is  shown.  To  obtain  the  bottom  block,  all  of 
the  top  blocks  must  be  removed.  Obviously,  if  it  is  desirable  from  a  system 
standpoint  to  have  the  bottom  block,  it  is  also  more  profitable  to  treat 
the  top  layer  as  ore  if  the  system  constraints  permit.  Using  the  static 
economic  cut-off,  the  top  level  would  always  go  as  a  waste. 


ore  waste 

ore  waste 

ore  waste 

-1  -3 

-1  -3 

-1  -3 

ore 

4 

FIGURE  1.4:  TYPICAL  CONDITION  OF  STATIC  ECONOMIC  CUT-OFF  FAILURE 

Also,  if  the  system  capacity  was  limited  to  only  one  more  block  and 
the  constraints  of  the  system  were  satisfied  by  any  one  of  the  four,  the 
most  profitable  selection  would  be  one  of  the  top  three. 


i 


14 


As  represented  In  [1],  there  has  been  some  attempt  to  incorporate  the 
system  influence  into  the  classifying  process,  but  none  that  considers  all 
the  influencing  factors. 

The  real  question  to  be  answered  in  cut-off  analysis  is  what  material 
to  mine  and  what  to  do  with  it  once  it  is  extracted,  so  as  to  maximize 
profits  subject  to  the  constraints  of  the  mining  system.  It  is  readily 
seen,  as  pointed  out  in  [14],  that  the  answer  to  this  problem  is  greatly 
influenced  by  the  state  of  the  mining  system,  which  varies  over  time.  Thus 
any  classification  decision  is  influenced  by  economies,  capabilities  of  the 
system,  assay  values  of  the  material,  mining  sequence,  desired  products,  and 
other  operational  constraints.  These  factors  may  interact  in  complex  ways 
over  time  and  hence  cause  a  variation  in  the  optimum  economic  or  grade  cut-off. 
Under  the  constraints  of  a  mining  system,  the  mining  plan  or  sequence  of 
mining  may  be  influenced  by  the  choice  of  cut-off,  and  hence  there  is  an 
interaction  between  these  two  elements. 

The  objections  mentioned  in  the  proceeding  paragraphs  are  eliminated 
by  the  linear  programming  model  to  be  discussed  in  the  following  chapters. 

This  model  provides  for  alternate  methods  of  classifying  blocks  and  considers 
the  system  state  in  its  decision  process.  Thus,  linear  programming  in  this 
sense  is  dynamic  marginal  analysis. 

To  obtain  a  better  understanding  of  why  the  model  provides  a  more 
optimum  cut-off  decision,  let  us  examine  the  cut-off  problem  from  the 
viewpoint  of  a  linear  programming  model: 

Maximize  cx  •=  z 
(1.4.1)  Subject  to  Ax  =  b 


x  >  0 


t 
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The  dual  to  this  problem  becomes; 

Minimize  jib  *>  v 
(1.4.2)  flA  >  c 

■n  unrestricted. 

The  linear  programming  approach  may  be  thought  of  as  decomposing 
a  system  into  a  number  of  elementary  components,  (activities),  describing 
the  interrelationships  of  the  activities  which  may  be  used  to  meet  the 
requirements  (b) ,  and  to  determine  the  level  of  activities  {x^}  ,  such 
that  the  requirements  are  satisfied  at  the  maximum  payoff.  In  this  sense, 
the  columns,  A^  ,  of  A  represent  the  input-output  coefficients  for 
activity  j  per  unit  of  activity  level.  For  example  a^  could  represent 
the  amount  of  requirement  b^  produced  by  activity  j  per  unit  of 
activity  level.  The  {c^}  can  be  considered  as  profits  per  level  of 
activity  J  . 

In  reference  to  the  primal  problem  (1.4.1)  the  dual  variables  arc 
called  multipliers  or  prices.  The  interpretation  of  the  dual  variables  as 
prices  comes  from  the  following  economic  interpretation  of  the  dual  problem. 

"Given  a  unit  profit  c^  for  each  activity  j  and  a  requirement  b^ 
for  each  resource  i  ,  what  must  be  the  unit  price  ,  of  each  resource 
such  that  the  total  value  of  the  resources  produced  or  consumed  by  J  is 
greater  than  or  equal  to  the  profit  c^  ,  and  such  that  total  value  of  the 
requirements,  irb  ,  is  minimal?"  [21]. 

The  condition  max  z  ■  min  v  is  a  consequence  of  the  duality  theorem, 
(6],  [21].  It  is  also  shown  in  [6]  that  if  nA^  <  c^  ,  it  is  profitable  to 
increase  the  level  of  activity  j  .  In  an  economic  sense,  what  is  occurring 
is  that  an  activity  whose  marginal  return  is  greater  than  its  marginal  cost 
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Is  being  substituted  for  an  activity  whom  marginal  re  turn  equals  marginal 
profit.  When  this  is  not  possible  with  present  prices,  the  optimal  solution 
has  been  reached. 


In  the  production  scheduling  model,  formulated  in  Chapter  2,  each 
activity  with  its  corresponding  level  X°?  represents  a  possible  method 
of  handling  block  (j)  in  the  mining  system  over  a  particular  time  period. 


Thus  we  have  a  set  of  alternatives  Qj  for  each  block  (j)  from  which  a 
convex  combination  of  alternatives  is  to  be  selected  and  combined  with 


other  convex  combinations  from  other  co  meet  the  requirements  in  such 

a  way  that  maximum  profits  arc  obtained. 

Using  a  static  cut-off  grade  or  economic  cut-off  selects  one  element 
from  Qj  without  considering  its  true  marginal  worth  to  the  system.  In 
selecting  an  element  of  in  this  manner,  one  is  never  quite  positive 

that  there  doesn't  exist,  another  element  e  such  that  nA^  < 

and  hence  profits  could  be  increased.  I.et  Z  be  the  maximum  profit  when 
the  are  restricted  to  one  element  and  Z  maximum  profit  with 

unrest  icted  Q^.  .  The  preceding  argument  shows  that  Z  <  Z  .  For  details 
on  linear  programming  the  unfamilar  reader  is  referred  to  [6]  or  [21]. 

It  is  now  clear  that  the  static  cut-offs  presently  used  in  mine  planning 
do  not  provide  for  maximum  profits  and  do  not  y ' e Id  the  best  possible  mining 
plan.  Yet  all  is  not  lost  since,  as  will  be  shown,  the  model  to  be  discussed 
in  this  paper  utilizes  a  dynamic  cut-off  which  fluctuates  with  the  system 
and  this  yields  a  mining  plan  which  is  as  close  as  presently  possible  to 


the  optimal  mining  plan  with  maximum  profits. 
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1.5  PREVIOUS  AND  KKLATKD  WORK 

The  use  and  potential  of  all  the  published  and  known  work  in  mine 
planning  has  been  limited  to  the  long  and  short  range  stages.  This  is, 
however,  closely  related  to  production  scheduling.  All  the  previous 
approaches  utilize  the  block  concept  except  the  work  by  Meyer  [16],  [17], 

The  static  cut-off  is  inherent  in  all  approaches  to  one  degree  or  another. 

Meyer  uses  a  pillar  approach.  He  formulates  the  problem  of  determining 
the  ultimate  pit  limits  subject  to  wall  slope  constraints  ns  a  mathematical 
programming  problem  with  a  nonlinear  objective  function.  Meyer's  formulation 
concentrates  on  the  geometric  problems  induced  by  the  slope  constraints. 

He  uses  the  principles  of  separable  programming  to  obtain  a  solution.  The 
pillar  approach  reduces  the  number  of  restrictions  and  variables  required, 
compared  to  the  block  concept,  and  hence  may  be  an  attractive  method  of 
determining  the  ultimate  pit  limits  when  accuracy  may  not  be  too  important. 

However,  the  pillar  approach  does  not  allow  for  any  extreme  variability, 
with  depth,  of  costs,  revenue,  or  geology,  and  thus  loses  accuracy.  When 
allowance  for  variability  is  introduced  the  formulation  converges  to  the 
block  concept. 

Meyer's  attempt  to  extend  the  pillar  concept  to  short  range  planning  is 
not  very  successful.  The  lack  of  variability  with  the  pillars  method  limits 
its  use  in  this  area.  His  proposal  is  to  maximize  profits  considering 
only  the  geometric  and  mining  capacity  constraints,  ignoring  completely  the 
flow  and  form  of  products  in  the  mining,  concentrating,  and  refining  process. 
Another  objective  is  the  prejudgement  of  what  is  ore  and  what  is  waste, 

(a  static  cut-off) . 


Mine  planning  at  the  Kcnnccotl  Copper  Corporation  a9  reported  in 
(8],  [19],  and  [20]  is  based  on  the  block  concept  with  simulated  extraction 
by  cones  and  is  patterned  after  the  traditional  methods.  Similar  to  most 
block  concept  systems,  Kennecott's  begins  with  a  preparation  of  a  mineralized 
block  inventory  and  an  economic  evaluation  of  the  blocks.  The  ore-waste 
evaluation  is  based  on  the  net  value  of  a  block  (a  static  economic  cut-off). 

The  Kennccott  open  pit  design  system  proceeds  toward  the  ultimate  pit 
limit  by  selecting  an  initial  trial  pit  (truncated  cone)  which  conforms  to 
the  desired  wall  slopes,  and  expanding  the  initial  cone  by  adding  conical 
increments  [20],  The  pit  limit  is  reached  when  no  additional  increment  can 
be  found  which  increases  the  total  net  value  of  the  material  with  an 
outline.  The  value  of  a  cone  or  increment  is  obtained  by  summing  the  values 
of  the  blocks  within  the  cones  or  increments.  One  of  their  special  problems 
is  the  variable  pit  slopes  that  must  be  maintained  in  some  of  their  operations. 
This  system  can  be  essentially  described  as  an  economic  evaluation  subject 
to  wall  slope  constraints. 

Kennecott’s  short  range  mine  planning  system  is  similar  to  its  long 
range  system,  as  it  is  patterned  after  the  traditional  manual  method  and 
subject  to  constraints  related  to  operating  slopes,  boundaries,  mine  system 
capacities,  and  operating  capabilities.  It  is  a  trial-and-error  method 
utilizing  a  trial  start  with  incremental  stages.  Each  increment  or 
combination  of  increments  represents  a  stage  in  the  development  of  the 
ore  body.  A  number  of  trial  plans  are  developed,  and  those  which  provide 
the  highest  dollar  value  with  desired  ore  volumes,  satisfactory  metal 
grades,  and  best  stripping  ratios,  are  selected.  This  method  does  not 
gnat  .mice  opt  i"  !itv,  a  1 1 hough  it  recognises  most  of  the  necessary  constraints 
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required  in  short  range  mining1,  planning  and  in  practical,  thus  gaining 
acceptance  by  the  mine  operators. 

Another  trial-and-error  block  concept  system  is  that  of  Minnesota  Ore 
Operations,  U.  S.  Steel  Corporation  as  reported  in  [1],  This  work  is 
probably  one  of  the  earliest  uses  of  the  block  concept  for  mine  planning. 

Here,  the  preparation  of  a  block  inventory  is  similar  to  what  has  already  bee 
described,  but  the  method  of  economic  evaluation  of  the  blocks  recognir.es  tic 
fallacy  of  a  pun  grade  or  economic  cut-off  bast’d  on  the  net  value  of  a 
block.  The  Minnesota  Ore  method  takes  cognizance  of  the  1  t  that  once  a 
block  of  material  is  removed  front  the  earth,  there  is  a  possibility  that  it 
is  better  to  treat  this  material  and  recover  sere  of  the  cost  rather  than 
completely  dispose  of  it  at  additional  cost.  The  objection  to  their  scheme 
arises  from  the  fact  that  this  classification  is  done  on  a  Mock  basis  and 
the  influence  of  the  entire  mining-refining  process  is  not  taken  into  account. 
This  method  also  considers  the  possibility  of  treating  an  ore  block  in  a 
number  of  ways.  (This  nay  not  be  possible  at  seme  operations.)  Thus,  in  a 
sense  this  approach  uses  a  dynamic  cut-off  on  a  single  bloc!  basis.  For 
mine  planning  and  scheduling,  a  number  of  blocks  are  ro-.ibin.  d  into  feasible 
mining  increments  called  shovel  units.  The  development  of  the  shovel  units 
is  accomplished  through  the  judgement  of  the  mining  engineer  with  full 
knowledge  of  the  material  classification  by  blocks,  equipment  requirements  of 
blocks,  grades  and  profit  values.  The  shovel  units  arc  the  entities  of 
scheduling  and  are  selected  on  a  profit  basis  and  combined  into  a  plan  or 
schedule  subject  to  mining  system  constraints.  Alternate  schedules  ire 
developed  to  meet  metallic  requirements  and  management  puli.  ie>;  mi  ,t  tri.il- 
and-error  bn  , i . 
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The  work  of  11.  Lorclis  and  I.  Grossman,  [15],  is  the  first  known 
nontrinl  and  error  approach  to  the  problem  of  doterming  the  ultimate  pit 
limits.  This  work  utilized  the  block  concept  and  their  assumptions  included 
equality  of  block  size  and  predetermination  of  what  was  ore  and  what  was 
waste.  The  Grossman-Lerchs  algorithm  is  a  directed  graph  algorithm  which 
divides  the  set  of  blocks  into  mineable  and  nonmincable  subsets,  positive 
and  nonpositive  subsets  on  the  basis  of  profitability  of  each  block.  Their 
algorithm  was  implemented  for  use  and  is  discussed  further  in  [10]. 

Although  their  primary  objective  was  to  determine  the  ultimate  pit 
contour,  they  did  recognize  that  there  are  many  ways  by  which  this  may  be 
reached.  To  accomplish  this  scheduling  they  suggested  the  use  of  an 
arbitrary  penalty  parameter,  by  which  the  profitability  of  each  block  could 
be  adjusted  to  determine  the  various  stages  of  development  towards  the 
ultimate  pit  plan.  This  penalty  parameter  could  be  thought  of  as  the 
assignment  of  Lagrange  multipliers  to  the  constraint  set. 

While  the  penalty  parameter  is  the  "germ"  of  a  good  Idea,  it  Is 
sometimes  a  difficult  task  to  arbitrarily  subset  the  proper  penalty, 
especially  when  numerous  constraints  are  present.  Also,  unless  some  method 
of  combining  plans  is  available,  the  method  does  not  always  converge  to 
the  optimum  plan  since  the  optimum  may  not  consist  of  entirely  whole  blocks. 

There  are  other  proposals  such  as  [7]  and  some  unpublished  work  at 
the  U.  S.  Bureau  of  Mines  which  has  contributed  to  the  area  of  mine  planning. 
However,  discussion  of  this  work  will  not  add  greatly  to  what  has  already 
been  given. 
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JL6_  TDK  CKO SSMAiN'-I.KRCHS  ALGORITHM 

KcrauHC  the  problem  to  which  the  Gros;  inan-I.erchs  algorithm  is 
Intended  to  be  applied  and  the  sub-problem  to  be  discussed  In  Chapter  3 
are  equivalent,  the  G-L  algorithm  will  be  presented  here  as  paraphrased  by 
Gilbert  [10]. 

1’ositivity  and  negativity  in  the  following  refers  to  positive  and 
negative  profits  of  a  subset  of  blocks. 

Algor i thm : 

Step  1.  Initialize:  each  block  of  the  whole  set  becomes  a  distinct  subset 
of  one  element. 

Step  2.  If  a  positive  subset  A  exists  which  is  constrained  by  a  negative 
subnet  B  ,  combine  both  to  form  a  new  subset  C  =■  AUB  .  1  f  no 

such  subset  A  exists,  go  to  Step  4,  otherwise  go  to  Step  3. 

Step  3.  Lxamine  only  subset  C  formed  in  Step  2  to  determine  if  an  advanta¬ 
geous  split  can  be  made.  A  split  will  be  advantageous  if  any  negative 
subset  D  exists  which  does  not  constraint  its  complement  in  C  . 

If  such  a  D  exists,  remove  it  from  C  and  return  to  Stop  2. 

Step  4.  Stop.  The  optimal  pit  is  identified  by  all  positive  subsets. 

For  the  original  presentation,  which  is  considerably  more  mathematical 
in  content  the  reader  is  referred  to  [15]. 

There  are  many  interesting  characteristics  of  this  algorithm  and  the 
problem  it  solves,  although  they  seemingly  were  not  recognized  by  I.erchs 
and  Grossman.  The  characteristics  of  the  problem  are  treated  extensively 
in  Chapter  3.  It  is  felt  that  the  method  to  be  proposed  in  Chapter  3  is  more 
efficient,  more  readily  implemented  and  more  adapt. <bl<  to  utvM'ul  blocks 
then  the  Grossnan-Lerchs  approach. 
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CilAlTKR  2 

FORMULATION  OF  GKNKKAI.  NATHIMATICAI.  MODKI. 

2.1  INTRODUCTION 

The  problem  to  be  i nvestignled  In  this  paper  ns  defined  in  Chapter  1  Is 
to  determine  an  "optimal"  production  schedule  for  an  open-pit  mining 
operation  for  a  planning  horizon  of  T  time  periods  (years,  months,  weeks, 
etc.).  The  term  "optimal"  refers  here  to  that  schedule  which  will  result  in 
the  maximum  total  profit  over  the  planning  horizon.  The  term  profit  may  be 
interpreted  to  bo  discounted  or  undiscounted  as  the  reader  desires. 

The  T  tine  periods  could  vary  in  length  depending  on  the  range  and 
purpose  of  the  scheduling  plan.  For  example,  one  could  view  the  problem  as: 
to  determine  a  yearly  schedule  for  the  next  ten  to  twenty  years  for  the 
purpose  of  future  economic  planning.  Another  possibility  is  to  find  a 
monthly  schedule  for  one  to  five  years  for  the  purpose  of  actual  production 
planning. 

As  will  be  shown  later,  the  techniques  to  be  presented  arc  also  appli¬ 
cable  to  the  problem  of  determining  the  so-called  optimal  pit  limits  used  in 
mine  evaluation.  Used  for  this  purpose,  the  method  would  provide  an  optimal 
ultimate  pit  limit  on  the  basis  of  present  technology  and  forecasted  economic 
conditions.  In  using  the  proposed  method  for  determining  the  long  range  pit 
limit,  consideration  can  be  given  to  factors  other  than  economic  (mining 
costs  and  revenue)  and  the  geometry  of  development,  such  as  the  form  and 
flow  of  products  in  the  mining-refining  system  and  other  factors,  which 
cannot  be  given  an  exact  numerical  value. 
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2.2  ASSUHn  IONS  OF  MODI il. 

2.2.1  Stat ciiii  nt  of  Assumptions 

In  order  to  obtain  a  workable  math cuat I cal  model,  a  number  ot 
assumptions  must  be  made  as  follows: 

1.  The  mineral  deposit  can  be  divided  into  a  finite  number  of 
mineable  units  called  blocks,  such  that  the  flow  and  form 
of  the  products  in  the  mining,  concentrating,  and  refining 
process  can  bo  described  in  terms  of  these  units. 

2.  All  blocks  are  the  same  size.  (Not  truly  essential  ,  but 
makes  development  of  procedures  easier.) 

3.  There  is  a  one-to-one  relationship  between  the  removal  of 
a  block  which  is  restricted  and  its  restricting  blocks. 

4.  It  is  possible  to  describe  the  allowable  mining  sequence, 
consistent  with  required  pit  geometry,  in  such  a  manner 
that  it  is  known  precisely  at  any  time  what  blocks  must  be 
removed  before  a  particular  block  may  be  removed.  Also 
this  sequencing  should  be  uniform  throughout  most  of  the 
probable  mining  area. 

5.  The  restrictions  for  mining  capacity,  milling,  capacity, 
required  volumes  of  products  at  certain  stages  of  the 
process,  required  grades  of  mill  feed  and  concentrate  and 
other  operating,  legal,  and  management  conditions  may  be 
expressed  as  linear  relationships  of  crude  volume. 

6.  One  is  able  to  assign  representative  values  to  each  block 
for  costs,  revenue,  volumes,  assays,  equipment  require¬ 
ments  and  other  data  necessary  for  the  relationships  in  ‘j. 


2 .2.2  Disc  uss  (.oil  of  Assumpt i ons 


Assumption  1: 

As  discussed  in  Section  1.3.2,  rectangular  blocks  provide  an  excellent 
means  of  describing  the  flow  and  form  of  products  in  the  almost  continuous 
mining,  concentrating  and  refining  process.  Actually,  open  pit  mining  Is 
essentially  accomplished  in  terms  of  blocks,  hence  this  is  not  a  very 
restrictive  assumption  for  our  model.  The  block  concept  has  been  used 
successfully  for  ore  estimation  in  a  number  of  cases  as  reported  in 
11,  8,  10,  15,  19,  20]. 

Assumption  2; 

Uniformity  of  block  size  yields  some  niceties  in  formulation,  however, 
it  is  not  a  necessity  for  the  treatment  of  the  problem  by  the  method  to  be 
proposed  as  will  be  seen  in  Chapter  3.  In  the  Interior  of  the  mineral 
deposit,  this  assumption  docs  not  limit  the  description  of  the  mining  process 
except  possibly  for  the  wall  s^opc  requirements.  This  exception  may  be 
overcome  by  varying  the  configuration  of  blocks  describing  the  allowable 
mining  sequence.  Another  objection  which  may  arise  is  that  it  may  be  more 
advantageous  for  planning  purposes  especially  where  selective  mining  is 
employed  to  vary  the  block  size  so  as  to  conform  more  closely  to  geological 
breaks.  However,  the  block  is  usually  quite  small  relative  to  the  entire 
amount  of  material  to  be  removed  in  a  planning  period  and  also  the  geological 
breaks  are  not  precisely  defined  (only  estimates  from  drill  hole  data). 
Therefore,  the  inclusion  or  exclusion  of  a  small  amount  of  undesirable 
material  in  a  block  will  not  be  detrimental  to  the  overall  optimal  schedule. 
It  would  seem  to  balance  out  in  the  final  plan. 

On  the  surface  boundary,  the  situation  is  somewhat  different.  Here  the 
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dime  its  ions  of  t  lie  I  op  Mucks  is  def  in.i  1 1  1  y  1 \  ■:  t  r  irled  by  1  In.  topography. 
For  c-x.im |»  1  •.  ,  .see  Figure  2.1. 
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FJCUKi:  2.1:  X-bl.CTION  SHOWING  r.lOCK  COI.I'HJIM'.ATJUN 

As  shown,  blocks  (1,1)  through  (1,R)  v;iry  in  heights  duo  to  the 
variation  in  the  surface  .  This  variation  can  easily  be  ovei  come  by  a  simple 
change  in  variable.  Let 


actual  volume  of  block  (i,j)  mined 
adjusted  block  volume  mined. 


Since  block  volumes  arc  assumed  equal  to  K  and  the  actual  is  if  we 

know  X,.  we  can  find  X,,  by  the  relation 
ij  ij 


x  -U 


Hence,  for  any  irregular  block,  the  above  transformat  ion  can  bo  used  and 
all  blocks  can  be  regarded  as  having  equal  size.  This  would  be  overwhelming 
if  we  had  to  do  it  on  all  blocks  but  with  only  surface  blocks  it  isn’t  too 
bad.  Of  course,  the  coefficients  of  the  various  constraints  which  include 
these  transformed  valuables  will  have  to  be  adjusted  also. 

A  reasonable  dimension  for  the  block  height  is  the  height  of  the  mining 
bench  of  the  operation  being  studied.  This  allow;  block  r»  ;< r* ••. > -n t - » t  ion  of 
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t ho  actual  Material  being  lui-n  d.  i.vn  i.  heights  in  present -  r.iinluj*  practice 
art*  in  the  range  of  40-50  feet. 

The  other  two  dimensions  of  a  bluet;  are  influenced  by  equipment 
capacities,  geological  structure,  wall  slopes,  manner  of  mining,  accuracy  of 
sample  data,  desired  use  of  block,  and  tin*  inconvenience  and  incapability  of 
manipulating  a  large  number  of  blocks.  It  is  geometrically  desirable,  due 
to  the  conical  shape  of  the  configuration  of  blocks  used  to  represent  the 
allowable  mining  sequence  to  make  the  other  two  dimensions  equal.  The  value 
of  this  dimension  will  be  called  w  in  all  further  discussion.  The  value 
of  w  can  usually  be  selected  within  the  following  limits: 

equipment  or  mining  method  bounds  w  <_  geological  bounds. 

Some  of  the  geometrical  conditions  on  w  will  be  discussed  in  relation 
to  the  block  configurations  representing  an  allowable  mining  sequence  since 
the  problems  are  related. 

Assumption  3: 

A  one-to-one  relationship  here  is  interpreted  to  mean  that  for  every 
unit  volume  which  is  mined  from  a  block  that  is  restricted  (overlaid)  by  a 
number  of  other  blocks  at  least  an  equal  volume  must  be  removed  from  each  of 
the  restricting  blocks.  For  example,  in  Figure  2.1,  assume  that  to  remove 
block  (2,4)  we  must  first  of  all  have  mined  blocks  (1,3),  (1,4),  and  (1,5). 
The  one-to-one  assumption  states  that  for  every  unit  volume  of  block  (2,4) 
mined  at  least  one  unit  must  be  removed  from  each  of  the  blocks  (1,3),  (1,4), 
(1,5).  It  is  conceivable,  therefore,  that  a  plan  would  result  in  only  taking 
1/4  of  each  block  (1,3),  (1,4),  and  (1,5),  and  1/8  of  bio. V  (2,4).  The 
question  will  arise  to  the  ’’operator”  how  to  get  1/8  of  lilo.k  (2,4).  If 
tiiese  blocks  occurred  on  a  boundary  of  a  mining  cut  other  than  the  surface, 
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it  i ii  entirely  possible  to  geL  such  answers  and  it  would  be  a  mineable  p'an. 
For  example,  a  fractional  block  plan  on  a  wall  slope  Jr.  shown  in  Figure  ?,2. 

Although  the  possibility  of  an  unuineable  situation  exists.  It  is 
thought  that  in  most  practical  problems  its  occurrence  will  be  negligible. 
The  rationale  for  such  an  opinion  is  based  upon  the  relatively  large  volumes 
that  will  be  removed,  the  occurrence  of  assay  value  trends  in  a  deposit,  and 
that  in  general,  all  else  remaining  relatively  fixed  costs  tend  to  increase 
with  the  depth  of  mining. 


Assumption  4 : 

The  allowable  mining  sequence  is  primarily  a  problem  of  gcomcLry.  In  an 
actual  mining  operation,  the  material  is  removed  in  such  a  way  that  the.  3a  fc 

wall  slope  is  never  exceeded  in  any  direction.  The  safe  wall  slope  is 

defined  as  the  angle  between  the  sides  of  a  mining  cut  ami  a  horizontal  plane 
at  which  the  material  will  stand  without  support  (angle  of  repos*.  1 .  The 
practice  of  benching  gives  the  open  pit  mine  a  slop-1  ike  structure.  (See 
Figure  1.1.)  These  conditions  justify  the  representation  of  the  outline  of  a 
mining  plan  or  the  volume  of  the  material  to  be  removed,  consistent  with  the 
wall  slope  constraints  and  mining  practices,  by  a  series  of  ’One-ffusLr’iv'i- 
cones  (see  Figure  2.3a).  Considering  the  vertical  sides  induced  by  the 

block  concept,  a  better  configuration  to  work  with  is  the  a'-sk-ccnc  shown  in 

Figure  2.3b. 

The  relationship  between  the  radii  of  the  respective  disks  is  given  by 


2 

r 


+  V?  +  r2 


whore 
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r  »  radius  of  disk  (cylinder)  which  lias  equal  volume  and  height 
as  f rustrum  of  regular  rone  with  radii  and  . 

Other  authors  dealing  with  related  problems  (10,  1T>)  have  proposed  a 
conical  representation,  but  when  the  block  concept  is  used,  the  representa¬ 
tions  arc  topologically  invarient.  That  is,  one  ends  up  using  the  disk-cone 
eventually  in  botli  approaches. 

At  this  point,  the  problem  is  to  determine  a  block  configuration  that 
closely  approximates  the  shape  and  volume  of  the  disk-cone.  A  desirable 
configuration  should  be  symmetrical  and  center  on  a  bottom  block.  The  block 
configuration  is  a  function  of  the  wall  slope,  0  ,  the  block  dimension, 
w  ,  and  the  height  h  .  Attempting  to  determine  a  block  configuration  and  a 
dimension  w  ,  which  satisfy  the  geometric  constraints  of  the  problem  and 
minimize  the  volume  difference  as  well  as  the  penalty  for  having  a  large 
number  of  blocks  gets  to  be  a  very  mossy  problem  and  it  is  not  the  function 
of  this  investigation  to  go  into  its  solution. 

As  we  have  bounds  on  w  (Assumption  2),  it  may  be  well  to  consider  w 
fixed  and  determine  a  block  configuration  that  satisfies  the  desired 
conditions.  This  can  be  accomplished  graphically  by  drawing  concentric 
circles,  representing  the  disks,  on  a  uniform  grid  of  w  by  w  squares  and 
comparing  the  volumes  generated  by  the  disks  and  various  patterns  of  blocks. 
Due  to  the  symmetry  and  uniform  grid  considerations,  the  patterns  will 
consist  of  combinations  of  a  9:1  configuration  and  a  5:1  configuration  [10J. 
These  configurations  are  shown  in  Figure  2.4. 

Two  of  the  possible  surface  patterns  generated  by  a  cubic  block 
configuration  scheme  of  1 : 5 : 9 : 5  for  an  overall  0=45  are  shown  in 
Figure  2.5.  The  configurations  and  volume.'.  generated  differ  due  to  the 
Starting  conditions  ( 1 : 9 : 5 : 9 : 5  or  1:5:9:5:9).  This  is  important  since  in  the 


v\  S  <  \\ 


RE  2.5: 


SURFACE  PATTERNS  FOR  1 : 5/1 : 9  CON'FICURAT  iON 


) 


33 


formulation  of  the  problem  we  will  want  to  assign  the  block  configuration 
representing  the  allowable  mining  sequence  a  priori  and  hence  may  enter  the 
configuration  at  any  level.  It  was  reported  in  [10]  and  confirmed  by  the 
author  that  a  1:5:9  configuration  gives  an  excellent  approximation  for  the 
volume  and  shape  of  a  diskical  cone  generated  by  a  regular  cone  with  a  side 

O 

slope  of  45  . 

The  allowable  mining  sequence  varies  from  operation  to  operation  due  to 
the  differences  in  permissable  wall  slopes  and  operating  practices.  Usually, 
the  wall  slope  is  constant  in  any  one  operation  but  there  are  cases  where  the 
wall  slope  varies  in  a  single  open  pit.  As  has  been  discussed,  the  allow¬ 
able  mining  sequence  for  any  open  pit  mining  operation  can  be  closely 
approximated  by  selecting  the  proper  configuration  and  block  size.  In  the 
cases  where  a  number  of  slopes  are  used  in  a  single  pit,  it  would  be 
advantageous  to  vary  the  configuration  rather  than  the  block  size.  In  the 
case  of  single  value  slopes,  it  appears  a  change  in  block  dimensions  would 

O 

be  the  best  thing  to  do.  For  example,  with  45  slopes,  w  ■  h  appears  to  be 
a  good  choice,  whereas  with  30  slopes,  w  ■  1.3  h  is  a  good  choice. 

Assumption  5: 

When  developing  a  model  for  open  pit  mine  production  scheduling, 
consideration  must  be  given  to  factors  other  than  the  geometric  conditions 
discussed  under  Assumption  3.  These  factors  provide  further  constraints  for 
the  model  and  may  be  classified  into  the  following  general  categories: 

1.  Form  of  the  material  within  a  block  and  the  method  of 
treatment  it  will  receive  in  the  mining,  concentrating, 
and  refining  system. 

2.  Equipment  availability  and  capacity. 
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3.  Avuilability  and  capacities  of  plants. 

4.  Manpower  uvnilubili ty. 

5.  Orderly  mining  of  the  pit, 

6.  Legal  and  physical  boundaries, 

7.  Desired  limits  on  assays  of  plant  feed. 

8.  Dcf.lred  form,  quantity,  and  assays  of  final  or  Intermediate 
products  of  the  system. 

As  will  be  illustrated  in  the  following  discussion,  most  of  the 
constralnts  in  the  above  categories  can  be  expressed  nr.  linear  relationships 
of  crude  volume. 

The  type  and  method  of  treatment  of  the  material  in  a  block  is  a  very 
Important  consideration.  All  too  often,  the  decision  as  to  whether  a  block 
ia  ore  or  waste,  and  If  either,  exactly  what  treatment  it  should  receive  in 
the  system,  is  made  without  considering  the  influence  of  the  entire  system 
and  other  blocks  in  the  system.  That  is,  the  material  is  ore  if  its  profit 
value  exceeds  a  certain  economic  cut-off,  and  Its  treatment  is  determined  on 
the  basis  of  metallurgical  tests,  which  are  valid  only  when  this  particular 
material  is  in  the  system  by  Itself.  A  priori  decisions  such  as  these  do 
not  always  lead  to  the  best  mining  schedule.  The  proposal  here  Is  to  let 
these  decisions  be  based  on  a  block’s  effect  on  the  profitability  of  the 
entire  system  operation  as  discussed  in  Chapter  1. 

This  proposal  can  be  included  in  the  model  by  letting; 


I 

r.p.t 


Yrpt 

*h 


*  amount  of  material  mined  from  block  h  , 
used  in  the  system  as  material  of  typo 
r  and  treated  by  method  (process)  ,  p  , 
during  period  t  . 


1.  *  total  volume  of  block. 

n 


35 


Consideration  of  equipment  availability  and  capacity  lends  to 
constraints  on  the  amount  of  material  that  can  be  taken  from  the  pit  or 
certain  areas  of  the  pit.  These  constraints  can  be  expressed  in  terms  of 
equipment  hours  per  ton  or  yard  of  material.  Equipment  used  in  an  open  pit 
mine  which  has  an  influence  on  the  mining  operation  capability  are  shovels, 
haulage  units,  drills,  and  other  auxillinry  mining  equipment.  Mining 
capacity  may  be  limited  in  a  certain  area  of  the  pit,  say  Q  ,  since  not  all 
the  equipment  can  be  concentrated  in  one  particular  area.  This  can  be 
expressed  for  a  single  period  as  follows: 

heQg 

where  *  hours  of  equipment  d  required  per  unit  volume  of  block  (h)  . 

Q  ,  g  ■  1,2,3,  ...  could  represent  a  number  of  small  areas  of  the  pit 

© 

or  it  could  be  the  entire  pit  operation,  in  which  case  the  above  constraint 
would  be  a  limit  on  the  total  capacity  for  equipment  type  d  .  Also,  the 
mining  capacity  may  be  limited  by  the  capacity  of  a  certain  sub-system  of  the 
mining  operation  to  handle  a  certain  amount  of  material  type  r  ,  for 
example,  r  ■  waste.  This  would  yield: 

1  1  X^pt  —  wa8te  handling  capacity  in  tons  or  yards. 

r*vaste  h 

If  an  operation  had  the  option  of  stockpiling,  there  would  be  similar 
constraints  for  this. 

The  capacity  of  plants  such  as  the  concentrator,  refinery,  and  even  the 
maintenance  facilities,  have  an  Influence  on  the  capability  of  the  entire 
system.  Examples  of  such  constraints  arc  as  follows: 


a^X^pt  <_  hours  available 


in  Q  for  each  type  of  equipment 

s 
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Manpower  is  a  very  important  element  in  n  ..'niiv,  opt-r  ■  The 

availability  of  r  anpov.er  n.ay  restriol  t  T  •  c.  p.v.  i  !  ; :  •  - ,  t  t !  or 

to  union  pressures  or  other  labor  relations,  un  i  for; :  i t  y  of  operation  may 
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have  to  be  considered.  There  are  a  number  of  operations  in  open  pit  mining 
which  can  be  expressed  as  functions  of  manpower  availability.  For  example, 
blasting  limitations  could  be  expressed  as: 


<  available  blasting  man-hours 


where  c  blasting  man-hours  required  per  unit  volume  of  block  (h)  . 
Although  not  directly  expressed  in  terms  of  man-hours,  lower  bounds  may  have 
to  be  considered  for  plant  and  equipment  availability  in  order  to  insure  a 
stable  environment  for  the  labor  force.  For  example. 


£  fjXT  _>  lower  bound  for  equipment  d  in  hours  in  period  t 

r,p,h 


where  f^  *=  hours  of  equipment  d  required  per  unit  volume  of  block  (h) 


£  b^X^1"  _>  lower  bound  for  sub-system  handling  material  r  wi 


th 


method  p  in  period  t 


where 

b^  =  hours  of  sub-system  per  unit  volume  of  block  (h)  treated 
in  process  p  as  material  r  . 

Conditions  other  than  the  geometric  factors  discussed  under  Assumption  3 
may  influence  the  orderly  mining  of  the  open  pit.  These  conditions  may  arise 
due  to  operating  procedures  such  as  access  for  transportation  units  or 
management  operating  policies.  Reasonable  expressions  for  constraints  of 
this  type  can  be  in  terms  of  ratios  of  material  on  different  levels  or  areas 
or  just  bounds  on  volume  from  particular  areas. 
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Legal  and  physical  boundaries  are  similar  to  the  orderly  mining 
constraints.  These  conditions  may  impose  limits  on  mining  volumes  and  on 
finAl  and  intermediate  products.  They  may  be  due  to  legal  contracts  such  as 
those  for  royalties  and  other  property  agreements.  Boundary  constraints  may 
also  arise  from  management  policies  and  locations  of  the  physical  plants. 

Desired  limits  on  assays  of  plant  feed  and  desired  form,  quantity,  and 
assays  of  final  or  intermediate  products  are  essentially  one  class.  They 
were  separated  in  listing  for  the  sake  of  emphasizing  the  importance  of 
control  of  plant  input.  Constraints  in  this  category  may  be  due  to  design 
characteristics  of  the  concentrator  and  refinery.  Also,  they  may  arise  from 
desired  products  as  influenced  by  marketability.  Design  characteristics 
give  rise  to  concentrator  feed  restrictions  expressed  as: 


Lower  bounds  on  feed 


h  i  i » 

p,h,reore 


where  q®  *  (crude  analysis-lower  assay  limit)  per  unit  volume  of  block  (h) 
for  element  s  ,  (s  could  be  iron,  silica,  copper,  etc.,  depending  on  a 
particular  operation),  or  lower  bound  on  the  assays  for  a  particular  treat¬ 
ment  of  type  r  material. 


T  n8PvrPt 

I  %  *h 

n 


1  0  . 


where  q^  ■  (crude  analysis-lower  assay  limit]  per  unit  volume  of  block  (h) 
used  as  type  r  material  in  process  p  . 

Upper  bounds  on  feed 

W" 

.  h  h 

p,h,rt.ore 


0 
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where  q8  -  (crude  analysis-upper  assay  limit)  per  unit  volume  of  block  (h)  . 
It  may  be  desirable  to  have  bounds  on  the  form  of  concentrator  input,  which 
can  be  expressed  as  ratios  of  different  materials. 


i  *jpt-  i 

P.h  p,h 


>  0 


where  vfs  expresses  the  desired  volumetric  ratio  of  materials  r  and  s 
Assay  bounds  on  the  concentrator  product  may  be  expressed  as: 

I 

p,h,r 


and 


l  i  ° 

p,h,r 


where 


n8**  *  [crude  assay-(upper  concentrator  assay  limit)  recovery] 
per  unit  volume  of  block  (h)  in  process  p  for 
element  s  . 

nSP  *=  [crude  assay-(lower  concentrator  assay  limit)  recovery] 
per  unit  volume  of  block  (h)  in  process  p  for 
element  s  . 


Similar  constraints  can  be  expressed  for  other  products  as  desired. 

It  is  the  opinion  of  the  author  that  the  constraint  areas  discussed  in 
the  preceding  paragraphs  cover  most  of  the  important  conditions  that 
influence  the  scheduling  plan.  Of  course,  there  may  exist  some  unusual 
conditions  which  are  not  adaptable  to  a  linear  expression,  but  it  is  felt 
that  they  would  be  of  minor  importance  in  determining  the  final  plan.  Many 


of  the  management  objection  which  arise  from  the  exclusion  of  such  unusual 
conditions  may  be  answered  by  sensitivity  analysis.  Part  of  the  discussion 
in  Chapter  3  will  cover  this  topic. 

One  condition  which  is  considered  important  by  the  operators  in  some 
cases  and  has  not  been  mentioned  in  any  discussion  so  far  is  minimum  pit 
bottom.  Some  operators  and  authors  have  expressed  the  opinion  that  the 
constraint  structure  must  provide  for  a  minimum  pit  bottom  area  in  order  to 
allow  wooing  space  for  equipment.  Even  though  it  can  be  implemented  quite 
easily  in  the  proposed  method  of  solution,  the  author  believes  that  it  is 
unnecessary  for  two  reasons:  (1)  because  of  geologic  trends,  it  is  felt 
relatively  few  blocks  in  a  plan  would  violate  this  constraint,  and  (2)  in 
actual  mining  practice,  the  expansion  to  the  minimum  full  working  bottom  is 
done  essentially  one  block  at  a  time.  Thus,  a  plan  which  bottoms  out  in 
only  one  block  can  be  thought  of  as  the  initial  stage  or  cut  into  the  lower¬ 
most  level. 

Assumption  6: 

Assignment  of  the  representative  values  to  each  block  required  for 
expressing  the  constraints  discussed  under  Assumption  5  has  been  discussed 
in  1.3.3.  As  stated  there,  a  number  of  assignment  procedures  have  been 
reported  in  [1,  8,  19,  20],  The  methods  are  dependent  on  the  particular 
operation  involved,  and  also  on  the  basic  data  available.  Great  care  should 
be  given  this  phase,  and  also  to  the  collection  of  the  basic  data  since  the 
solution  of  the  scheduling  problem  is  only  as  good  as  the  data  on  which  it 


is  based. 


2.3  MATHEMATICAL  FORMULATION 


Based  on  the  preceding  assumptions,  the  problem  of  determining  a 
extraction  schedule  for  an  open  pit  mine  over  a  duration  of  T  time 
which  will  yield  the  maximum  profit  may  be  formulated  as  a  large  srd 
linear  programming  problem  with  the  following  structure; 


Maximize 

cix1 

+  c2X2 

T  T 

+  .  .  .  +  c  X 

Subject  to 

D1X1 

+  d2x2 

T  T 

+  ...  +  D  X 

AjX1 

A,X2 
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Xfc  *  X^pt  RPH  x  *  component  vector  for  each  period  t  ■  1,  ...,  T  . 

R  ■  number  of  material  types. 

P  ■  number  of  possible  methods  of  treatment. 

H  ■  IJK  »  number  of  blocks. 

I  ■  maximum  level  dimension. 

J  »  maximum  column  dimension. 

K  =  maximum  section  dimension. 

Aj  »  Bi  x  n  matrix  corresponding  to  the  coefficients  of  the 

constraints  discussed  under  Assumption  5  for  period  1  .  n  ■  RPIJK 

bi  "  mi  x  1  vector  correspond  <g  to  the  right-hand  side  of  these 
constraints. 

0  -  w  x  1  zero  vector 

«  mQ  x  n  matrix  corresponding  to  the  pe  lod  i  portion  of 
coefficients  of  constraints  discussed  under  Assumption  5  which 
relate  to  the  entire  planning  duration. 

d  ■  mQ  x  1  vector  corresponding  to  the  right-hand  side  of  these  con¬ 
straints. 

c*1  «  1  x  n  vector  of  cost  or  profit  coefficients  for  period  p  . 

■  volume  of  block  (h)  . 

Bit  ■  w  x  n  matrix  of  coefficients  for  the  constraints  which  describe 
the  allowable  mining  sequence  as  discussed  under  Assumption  4. 

For  the  purpose  of  explaining  the  detailed  structure  of  the  matrices 

it  is  convenient  to  change  from  the  single  subscript  notation 
h  ■  1,2,  ...,  H  for  $ach  block,  to  a  more  natural  3-dimensional  notation 
i  *  1,2,  ...,  X  ,  J  *  1,2,  ...,  J  and  k  ■  1,2,  ..., 


K  ,  as  shown  in 


levels 


i 


columns 


J 


FIGURE  2.6:  BLOCK  STRUCTURE 


now  are  designated  by 


Figure  2.6.  Our  block  volume  variables 

Yrpt 
ijk  ’ 

Consider  t lie  sequence  extraction  constraints  for  block  (232)  (directly 
below  block  (132))  assuming  we  must  remove  the  five  blocks  restricting  it 
as  discussed  under  Assumption  4.  Also  assume  the  schedule  duration  to  be 
two  periods,  blocks  to  be  classifiable  as  ore  or  waste,  and  that  it  is 
possible  to  treat  the  ore  by  two  different  methods.  Then  we  have,  in 
accordance  with  Assumptions  3,  4,  and  5,  the  following  constraints  which 
illustrate  the  structural  form  of  the  B.  matrices: 
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PERIOD  1  PERIOD  2 


These  relationships  express  the  fact  that  to  get  at  block  (232)  in 

Period  1,  blocks  (122),  (131),  (132),  (133),  and  (142)  must  be  removed  in 

Period  1  irrespective  of  what  is  done  with  the  material  after  it  is  mined. 

To  obtain  block  (232)  in  Period  2,  the  sequence  of  blocks  (122)  -  (142) 

can  be  removed  in  Period  1  or  Period  2. 

It  is  easily  seen  from  the  coefficient  matrix  of  (2.3.2)  that  each 

B,  can  be  expressed  as: 

ip 

B .  =  [EE  E  E] 

ip 


by  re-arrangement  of  the  columns,  where  each  E  is  a  matrix  with  exactly 
one  -  1  and  one  +1  in  each  row,  and  all  other  elements  are  zero. 

Thus  E  is  the  transpoce  of  the  node  arc  incidence  matrix  of  a  network- 
flow  problem,  [see  Appendix  or  [2],  [21]] 

Before  giving  a  theorem  which  will  demonstrate  an  important  advantage 
of  this  structure,  consider  the  constraints: 


I  l 

r,p  r,p 


+  l 

r.P 
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ijk  ' 


Dividing  each  of  these  constraints  by 


luk  gtves 


l 

r»P 


I 

r*P 


+ 


l 

r,P 


1 


where 
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This  same  transformation  can  be  carried  out  for  each  constraint  of 
2.3.1  and  is  possible  because  of  Assumption  2,  (t^^  *  K)  .  Note  the 

matrices  ,  and  remain  unchanged,  but  the  b^  and  d  are 

bi  -  d 

transformed  to  b.  =  - -  and  d  *  - -  .  In  all  further  discussion, 

1  ijk  ijk 

problem  (2.3.1)  will  be  considered  in  this  transformed  form  unless  stated 
otherwise. 


Theorem  2.1: 


If  the  Dantzig-Wolfe  Decomposition  technique  is  used  to  solve  (2.3.1) 


with 


B 2 /  +  B22x2 


(2.3.3)  Bt1X1  +  Bt2X2  +  . 


. .  +  bttx  <  0 


l  XhPl  +  l  XhP2  +  ‘ 

r  »P  r,p 


•  •+  I  <pT :  i 

r,P 


for  all  ijk 


x*pt  >  0 

n  - 


as  the  sub-problem  constraintSj  then  the  actual  sub-problem  to  be  solved 


can  be  reduced  to: 


-11  -2  2  -T  T 

Maximize  d  K.  +  d  Y,  +  . .  .  +  d  Y, 
h  h  h 


(2.3.4) 


Subject  to  EY 


EY1  +  EY2 


12  T 

EY  +  EY  +  ...  +  EY  <  0 


where 


12  T 

Y.+  Y.  +  . . .  +  Y*  <  1 
n  h  h  - 


yJ  >  o 

n  - 


t  •  C7 


-sut  -rpt 

c,  *>  max  c,  r 
h  h 
r  *P 


nt  -sut  sut  „hsut  .hsut 

d,  *=  c  «=  c,  ••  QD  -  tt  A 
n  n  n  t  t  t 


Proof ; 


From  the  structure  of  each  *=  {E,E . E}  ,  it  is  seen  that  only 

one  of  the  rp  variables  X^*1  can  appear  in  the  basis  solution  of  (2.3.3), 

otherwise  the  basic  columns  would  be  dependent,  contradicting  the  Independence 

su  t 

of  a  basis.  If  some  X^  does  occur  in  the  basic  optimal  solution  at  a 

positive  level,  it  must  correspond  to  c®ut  *  max  c^C  .  This  is  true  since 

h  h 

r  >P 

the  columns  corresponding  to  the  set  of  variables  X^*,r  ■  1,  ...,  R  , 
p  ■  1,  ....  p  ,  are  all  the  same  and  at  optimality 


n  *sut 
0  >  c, 

“  n 


T 

-sut  r  -  Dhsut  ^  -rpt  c  -  _hrp 

c,  -  )  u.B_  >  c,  K  -  )  tt.B„  v 

h  l  It  -  h  6  l  It 


for  all  r,p  . 


» 
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Hence,  to  solve  problem  of  (2.3.3),  we  need  only  consider  those  columns 
corresponding  to: 

vt  sut  /  -sut  -rDt 

Yh  "  Xh  /  ch  "  max  c< 

ru' 

Doing  this  we  get  (2.3.4). 

Q.E.D. 

Theorem  2.1  will  be  used  to  great  advantage  In  the  chapters  to  follow 

since  It  provides  considerable  simplification  to  problems  In  the  form  of 
(2.3.3). 
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CHAPTER  3 

THE  SINGLE  PERIOD  PROBLEM 


3.1  INTRODUCTION 

In  order  to  investigate  the  properties  of  the  scheduling  problem  as 
formulated  in  Chapter  2  and  to  simplify  the  development  of  a  solution 
technique  for  a  T-period  problem,  the  single  period  problem  will  first  be 
considered. 

The  one-period  problem  has  the  following  form: 


(3.1.1) 


Maximize 
Subject  to 


cX 


AX 

BX 


l 

r.P 


*  b 

<  0 

V 

<  1 

E 


:°J 


for  all  ijk 


The  single-period  problem  is  similar  to  the  ultimate  pit  limit  problem 
of  which  a  number  of  reports  have  appeared  in  the  literature  as  Indicated  in 
Chapter  1;  [1,  5,  10,  16,  17,  18].  The  problem  as  considered  by  these 
authors  excludes  the  constraint  set  AX  ■  b  and  considers  only  an  apriori 
classification  of  the  material  which  is  essentially  setting  B  *  E  .  (E  ■ 
transpose  of  node  arc  Incidence  matrix,  see  Appendix).  This  type  of  an 
approach  yields  a  rough  idea  of  the  mineable  reserves  and  is  valuable  in  a 
preliminary  evaluation  study  of  a  deposit.  However,  it  is  the  opinion  of 
the  author  that  it  is  more  desirable  to  study  the  "reserve  picture"  under 
the  influence  of  the  constraints  which  limit  the  form  and  flow  of  products 
in  the  entire  system. 

The  approach  presented  in  this  paper  is  applicable  to  the  ultimate  pit 


* 


b) 


limit  problem  as  well  as  many  extentions  which  may  be  valuable  in  the 
evaluation  and  design  of  a  mining,  concentrating  and  refining  system.  For 
example,  it  may  be  important  to  know  the  ultimate  pit  limits,  or  reserve 
tonnage  for  a  particular  time  horizon  when  volume  restrictions  are  imposed  on 
various  components  of  the  system  and/or  when  the  average  assay  of  mill  feed 
or  concentrates  are  limited  to  a  particular  range.  Markets  for  the  system 
products  may  also  impose  constraints  which  will  influence  the  total  reserves. 
There  are  many  other  problems  which  may  be  investigated  by  the  methods  of 
this  paper  but  those  presented  seem  sufficient  to  point  out  the  versatility 
of  the  approach. 

The  linear  programming  problem  (3.1.1)  will  in  most  practical 
situations,  exceed  the  capabilities  of  presently  available  linear  programming 
codes  and  computing  facilities. due  to  the  dimension  of  the  B  matrix. 

(Recall  from  the  discussion  in  Chapter  2  that  there  are  approximately  seven 
constraints  in  BX  *  b  per  block.  The  A  constraints  will  be  relatively 
few  in  number  (less  than  40).)  Because  of  these  limitations,  it  seems 
desirable,  even  necessary,  to  seek  other,  more  compact,  methods  of  solution. 


52 


3.2  ANALYSIS  OF  THE  PROBLEM 
3.2.1  Decomposition  of  Problem 

Due  to  the  etrueture  of  problem  (3.1.1),  the  relatively  small  else  of 
the  A  matrix  and  the  simplified  structure  of  the  B  matrix  which  can  be 
reduced  by  virtue  of  Theorem  2.1,  it  appears  attractive  to  examine  the 
application  of  the  Dantrig-Wolfe  decomposition  principle  [6]  to  this  problem 
Applying  the  decomposition  principle  to  problem  (3.1.1)  using 

the  constraint  set 

master  problem: 


BX 


<  b 


l  *15.  :  i 


r.p 


ijk 


as 


crp 

ijk 


>  0 


of  the  subproblem  yields  the  following  full 


(3.2.1.1) 


Maximize 

Subject  to 


where  JT  is  an  extreme  point  solution  of  the  subproblem.  For  optimality 
it  is  required  that  cX*  -  wAX*  -  s  <  0  for  all  i  .  Here  ir  is  the 

L  i 

1  x  m  vector  of  multipliers  for  the  equations  £  A  .AX  •  b  and  s  is 

i-1  1 

the  scalar  multiplier  corresponding  to  the  equation 
L 

l  A  -  1  [6]. 

i-1  1 


Therefore,  the  subproblem  becomes 
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:  • 


t 


Maximum 
Subject  to 


(3.2. 1.2) 


(c  -  nA)X 
BX  <  0 

m 

y  xrP  <  i 

r  D  2 

r»p 


V  ljk 


xrP 

*ijk 


>  0 


The  general  procedure  (aa  given  in  (5,  6])  assuming  we  have  at  hand 
an  Initial  feasible  solution  to  the  restricted  master  (3. 2. 1.3) 


Maximum  l  X^cX' 
i-1 


(3.2. 1.3) 


AX  ■  b 

-  1 


X.  >0 

1  * 


i;i . . 


is  to  determine  it  ,  s  by  solving 


wAX1  +  s  -  cX1 


for  1  ■  1,  . . . ,  q 


Then  solve  the  subproblem  (3. 2. 1.2)  to  determine  if  maximum  cx  -  nAX  <  s  . 
If  this  holds,  then  the  optimal  solution  to  problem  (3. 2.1.1)  has  been 
obtained  and  hence  to  problem  (3.1.1)  through 


f  v 

I  * 


x  -  l  \X 

i-i 


If  maximum  cX  -  itAX  >  s  then  it  is  profitable  to  Introduce  the  new  vector, 
^AXj  with  cost  coefficient  cX  ,  into  the  basis. 

After  introducing  the  new  vector,  new  it  ,  s  are  determined  and  the 
process  is  repeated. 


i 


f 


«e 
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Once  the  subproblem  has  been  solved,  introducing  a  new  vector  into  the 
master  problem  and  solving  it  is  relatively  simple,  since  the  master  is 
rather  a  small  (less  than  40  constraints)  linear  programming  problem. 


3.2.2  The  Subproblem 

The  problem  remains  how  to  solve  the  subproblem  (3. 2. 1.2)  since  as  it 
appears  it  is  a  large  linear  programming  problem.  Fortunately,  it  can  be 
reduced  by  employing  Theorem  2.1.  By  virtue  of  this  theorem  the  new  sub* 
problem  becomes 


Maximize 

(c  -  xA)*Y  -  cY  -  Z 

(a) 

(3.2. 2.1) 

Subject  to 

EY  <  0 

m 

(b) 

0<  Y<  1 

(O 

where  c  »  (c  -  nA)*  ■  (c8“k  "  *Ai^8UJ  such  that 


su  .ljksu 
-  wA  j 


cijk  ‘  cyk  *  ’* --  ‘  (cIjk  - 

r*P 


irA 


ijkrp^ 


,su 


and  Y  -  {YiJk}  such  that  Y^  -  «1Jfc 


which  corresponds  to  c 


IJk  ’ 


Lemma  3. 2. 2.1: 

There  always  exists  an  optimal  solution  to  the  problem  (3. 2. 2.1). 

Proof : 

A  feasible  solution  Y  ■  0  always  exists  and  any  feasible  solution  is 
bounded  (Z  ^  ®)  .  Therefore  as  a  consequence  of  the  Duality  Theorems  of 
linear  programming  (see  (5,  21]),  it  is  known  that  the  dual  of  this  problem 
always  has  a  feasible  solution,  and  an  optimal  solution  to  both  problems 
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exists  such  that  max  Z  ■  min  (the  value  of  the  dual  problem). 

Q.E.D. 

t 

Due  to  Lemma  3. 2. 2.1  there  need  be  no  concern  for  any  X^  >  0  other 

than  those  which  satisfy  £  X.  *  1  in  the  master  problem. 

1  1 

The  dual  to  problem  (3. 2. 2.1)  Is: 

Minimize  P  -  l  p...  (a) 

ijk 

Subject  to  uE  +  pi  >  c  (b) 

u,p  >  0  (c) 

It  is  Immediately  clear  that  a  feasible  solution  to  problem  (3. 2.2.2) 
la: 


u  -  OtPjjk  "  ma*  (°*cijk)  for  a11  * 


Lemma  3. 2. 2. 2: 


The  following  bounds  exist  for  max  Z  of  problem  (3. 2. 2. 1)  : 


II  0  (a) 

:  l  \  <»> 

h/v° 


Proof; 

Part  (a)  Is  obvious  since  max  Z  >  cY  for  any  feasible  Y  and  Y  •  0 
Is  feasible  for  (3. 2. 2.1)  hence  max  Z  >  0  . 

To  show  (b) ,  the  weak  duality  theorem,  [5,  21],  states  that: 

P  ■  l  >  cY  *  Z  for  any  feasible  p  ,  Y  ,  and  u  .  Hence  setting: 
h 


u  -  0 

ph  -  0  lf  'hi° 


and 


ph  ■  ch  lf  ch  *  0 

results  In 


K  -  l  s  i  'Y 

h/ch>0 

for  all  feasible  Y  and  hence  for  an  optimal  Y  corresponding  to  max  Z  . 

Q.E.D. 


The  necessary  and  sufficient  optimality  conditions  for  problems 
(3.2. 2.1)  and  (3. 2. 2. 2)  are: 


EY  <  0 

m 

0  <  Y  <  1 

■  ■  / 

_  i 

uE  +  pi  >  c 
u  >  0  ,  p  >  0 
u(EY)  -  0 
p(I  -  Y)  -  0 
(c  -  uE  +  pI}Y  »  0 


Primal  Feasibility 


Dual  Feasibility 


Complimentary  Slackness 


3.2.3  Properties  of  Primal  and  Dual  Suboroblems 

An  Investigation  of  the  properties  of  the  subproblem  will  lead  to  a 
better  understanding  of  the  problem  and  aid  In  stating  an  efficient  method 


of  solution. 


5? 


As  was  pointed  out  in  Chapter  2,  the  structure  of  the  matrix  E  in 
problems  (3. 2. 2.1)  and  (3.2.2. 2)  has  the  characteristics  of  the  transpose  of 

'•v 

a  node-arc  incidence  matrix  of  a  network  or  graph, at  most  two-*nonzero 
elements  per  row. 

Examining  the  dual  problem  (3. 2. 2. 2),  it  is  seen  that  each  block  (ijk) 
corresponds  to  a  node  in  the  network.  The  arcs  corresponding  to  the 
variables,  p  ,  are  incident  out  from  the  set  Y  -  {(ijk)}  .  The  arcs 
corresponding  to  the  variables,  u  ,  are: 

Incident  into  node  (ijk)  if  the  coefficient  of 
u  in  the  E  matrix  is  -1 

q 


and 


Incident  out  from  node  (ijk)  if  the  coefficient  of 
u  in  the  E  is  +1  . 

q 

The  variables  p  and  u  can  be  thought  of  as  flow  variables  on  their 
respective  arcs.  Taking  the  dual  constraint  set  and  adding  nonnegative 
slack  variables  yields: 

uE  +  pi  -  vl  ■  c  , 


the  conservation  equations  of  the  network  corresponding  to  the  dual  problem. 
The  variables  v  *  {v^^  correspond  to  flow  on  arcs  leading  into  each  node. 
Each  component  c^  of  c  corresponds  to  a  given  quantity  of  flow  into  or 
out  of  each  node  depending  on  the  sign  of  c^jjc  • 

This  assignment  is  as  follows: 


I 
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If 

If 


cuk  * 0 

cljk  i  0 


*0 

■> 


flow  Into  node  (ijlc)  -  a  demand  of 

flow  out  of  node  (ijk)  -  a  supply  of  . 


To  clarify  these  relations,  consider  the  small  two  dimensional  system 
of  blocks  shown  in  Figure  3.1.  The  block  number  appears  In  the  upper  left 
corner  and  c^  at  the  bottom  of  each  block.  Assume  the  allowable  mining 
sequence  is  such  that  for  any  block  (lj)  or  partial  block  removed  from  a 
level  (other  than  level  0)  the  block  directly  above  it  (1-1, j)  and  its 
neighboring  blocks  (1-1, j-1)  and  (i-l,j+l)  must  also  be  removed. 

The  network  corresponding  to  the  dual  problem  for  this  system  of  blocks 
is  shown  in  Figure  3.2. 

The  first  thing  one  notes  when  looking  at  the  network  is  that  there  are 
no  arcs  corresponding  to  the  variables  {p^^l^O)  and  {v^.J-3}  .  In 
other  words,  these  variables  have  been  set  equal  to  zero  and  their  corres¬ 
ponding  arcs  Ignored.  This  is  a  general  characteristic  of  the  networks 
represented  by  the  dual  problem  and  is  the  result  of  the  following  lemma. 


Lemma  3. 2. 3.1: 
The  problem 


Minimize  J  p 
ijk 

Subject  to  uE  +  pi  -  vl  ■  c 
u,p,v  >  0 

has  an  optimal  solution  that 


o 

■ 

for 

all 

pijk 

-  0 

for 

all 

v. 

ijk 

ijk 

i  i  0  and 

i  i  max  i  if  there  exists 


t 


S^mmxgsk 

SMSas9M9 


FIGURE  3.2:  NETWORK  REPRESENTATION  OF  BLOCK  SYSTEM  IN  FIGURE  3.1 
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a  "free  path "  from  node  ( ij k)  to  node  (i-ltjtk).  A  " free  path"  ie 
defined  ae  a  path  for  which  there  ie  no  associated  cost. 

Proof : 


From  Lemma  3. 2. 2.1,  it  Is  known  that  the  above  problem  always  has  an 
optimal  solution.  Assume  that  there  exists  and  optimal  solution  with 
P  -  {iJk/pijk  >  0,1*0}  and  P°  -  {Ojk}  and  V+  -  {ljk/vijk  >  1  *  max  1} 
and  Vn  ■  {max  l,j,k}  . 


From  the  statement  of  the  lemma,  all  nodes  are  accessable  from  the  set 

Vn  by  a  free  path.  Hence  for  some  par  e  V+  so  that  v  ■  e  .  Set 

pqr 

v  *0  and  for  some  (stv)  c  Vn  such  that  (stv)  >  (pqr)  c  V+  .  Set 
pqr 

vstu  "  ^stv  +  e  and  chan6e  the  alonB  the  path  from  (8tv)  to  (pqr) 
by  +e  In  each  arc.  Also  It  is  clear  that  there  exists  a  free  path  from 
all  (ljk)  c  P+  to  P®  .  Hence  setting  p^k  -  0  for  (ijk)  t  P+  and 

p0pr  "  fiQpr  +  cijk  for  >  ^0pr^  and  lncrcasln8  flow  by  along 

this  path  gives  the  desired  result. 

Q.E.D. 


It  should  be  noted  that  this  problem  could  be  treated  as  a  minimum  cost 
flow  problem  and  solved  by  use  of  the  out-of-kilter  algorithm.  This  does 
not  seem  particularly  attractive  with  the  present  form  of  the  problem. 


^Accessability:  If  n^  and  n^  are  nodes  such  that  there  is  a  path  from 

n^  to  n^  then  n^  is  "accessable"  from  n^  written  as  n^  >  n^  .  The 

accessable  set  D(n^)  is  the  set  of  all  nodes  accessable  from  n^  , 

(Oj  c  D(n^)  .  Define  I)(n^)  *  D(n^)  +  n^  .  The  work  "path"  used  in  the 

above  explanation  is:  a  sequence  of  arcs  {a.,  . ...  a  ,  a  ,,  ...,  a  } 

1  1  1+i  n 

such  that  the  terminal  node  of  each  arc  coincides  with  the  initial  node  of 


the  succeeding  arc;  i.e.,  all  arcs  in  the  path  have  a  similar  orientation. 
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however,  due  to  the  large  number  of. arcs  one  would  have  to  consider. 

The  primal  variables,  are  called  node  numbers  or  node  potentials 

in  their  relationship  to  the  network  formulation  of  the  dual  problem. 

Yljk  "  *  mRans  a  block  is  to  be  mined  while  Y^  ■  0  means  the  block 
(ljk)  will  be  left.  As  will  be  shown  later  (Theorem  3.3),  the  solution  to 
the  subproblem  is  in  terms  of  blocks  to  be  mined  (Y^^  ■  1)  or  left 
■  0)  .  It  will  be  clear  that  the  subproblem  provides  a  system  of 
blocks  to  the  master  problem  which  selects  and  combines  the  systems  into  the 
optimal  mining  plan. 

The  theorems  to  follow  characterize  the  solutions  to  the  primal  and  dual 
subproblems,  (3. 2. 2.1)  and  (3. 2. 2. 2)  and  hence  aid  in  presenting  a  method  of 
solution. 


Theorem  3.1: 

If  for  any  node  corresponding  to  a  block  which  is  not  restricted  by 
any  other  block ,  the  condition  o^  >  0  holds,  then  p^>  0  and  1^  »  l  in 
the  optimal  solution  of  problem  3.2. 2.1. 


Proof: 


Looking  at  the  dual  equations  for  any  node,  h  ,  which  is  not 

restricted  by  other  nodes  (such  as  the  nodes  on  level  "0"  in  Figure  3.2) 

it  is  seen  that  the  arcs  incident  on  h  ,  excluding  the  arcs  corresponding 
to  p^  and  ch  are  all  incident  into  node  h  .  Mow  if  c^  >  0  (incident 

into  h)  the  p^  is  the  only  possible  flow  leaving  h  .  Kence  from  the 

dual  inequality  for  node  h  , 


u<i.  +  Pk  2 
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it  is  clear  that  p^  >  c^  >  0  .  From  the  complimentary  slackness  condition 
p^(l  -  Y^)  *=  0  it  is  clear  that  p^  >  0  =>  Y^  *=  1  . 

Q.  E.D. 

Theorem  3.1  provides  a  means  of  classifying  unrestricted  blocks  by 
immediate  inspection  of  ch  .  There  may  be  blocks  in  level  "0"  that  are 
to  be  selected  (Y^  *  1)  for  other  reasons  which  will  be  brought  out  by  the 
theorems  to  follow. 

Theorem  3.2: 

If  some  ,  then  every  corresponding  to  (k)  t  D(h)  is  l 

in  the  optimal  solution . 

Proof : 

Let  D.  (p)  ■  tq/(p,q)  e  A  >  .  Where  A+  ■  set  of  arcs  incident  out 

1  l  Pi  p 

from  node  p  .  From  the  primal  equations:  EY  <  0  we  have  that  if 

»  1  then  Y^  ■  1  for  k  e  D^(h)  .  Then  consider  i  c  D^[D^(h)]>  again 
from  the  primal  relationships  Y^  ■  1  .  Continue  in  this  manner  until  the 
nodes  on  level  "0"  are  reached  for  which  the  accessable  set  is  empty  and 
the  theorem  is  proven. 

Q.E.D, 

The  next  theorem  is  important  in  that  it  shows  that  we  need  only  concern 
ourselves  with  potentials  {Y^^}  0  or  1  .  Thus  the  primal 

subproblem  solutions  are  sets  of  zeros  and  ones  as  indicated  previously. 

Before  giving  the  next  important  theorem,  we  must  discuss  briefly  the 
unimodular  property  of  matrices.  A  matrix  A  is  called  unimodular  if  the 
determinant  of  every  square  submatrix  of  A  equals  0  ,  +1  or  -1  . 

The  following  lemma  is  given  by  Hoffman  and  Kuhn  in  [13]. 


Lemma  (Hoffman  and  Kuhn): 


If  in  a  system  of  linear  inequalities  with  integral  coefficients  and 
constant  terms,  every  nonsingular  submatrix  of  the  coefficient  matrix  has 
determinant  +1  ,  then  every  extreme  solution  is  integral. 

T 

Since  It  is  well  known  that  every  matrix  [E  I]  la  unlmodular  If  E 
is, the  next  theorem  follows  Immediately. 

Theorem  3.3: 

«  l  or  0  in  the  optimal  solution  of  the  subproblem: 

Maximize  cY 
Subject  to  EY  <  0 
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3.3  EQUIVALENT  FORMULATIONS  OF  SUBPROBLKM 

The  subproblems,  (3. 2. 2.1)  and  (3. 2. 2. 2),  analyzed  in  Section  3.2  were 
stated  as  follows: 


Primal:  Maximize  cY  ■  Z  (a) 

(3.3.1)  Subject  to  EY  <  0  (b) 

0  <  Y  <  1  (c) 

H 

Dual:  Minimize  P  ■  \  p.  (a) 

h=l  n 

(3.3.2)  Subject  to  uE  +  pi  >  c  (b) 

u,p  >  0  (c) 


From  a  network-flow  point  of  view  the  objective  in  solving  the  dual 

problem  (3.3.2)  is  to  induce  flow  into  the  network  in  such  a  way  that  the 

requirements  of  the  c^  are  satisfied  and  as  little  as  possible  goes  through 

the  arcs  corresponding  to  the  p^  variables.  As  indicated  in  Section  3.2.3, 

flow  is  forced  into  the  network  by  the  positive  c^  and  taken  from  the 

system  by  the  negative  c^  .  For  any  node  (h)  there  are  only  two  ways  for 

flow  that  is  forced  in  (condition  when  c^  +  \  u^  *  0)  to  get  out, 

i/heD1(i) 

one  way  is  to  route  flow  along  an  arc  corresponding  to  u..  ,  l.e.t  to 

nk 

another  node  k  ,  or  remove  flow  via  a  p^  arc.  Therefore,  flow  leaves  the 
system  (or  reaches  the  sink)  only  through  the  p^  arcs  or  through  the  arcs 
corresponding  to  c^  <  0  .  This  indicates  we  may  replace  all  paths  from  a 
node  (h)  such  that  c^  >  0  to  a  node  (k)  v/ith  <  0  by  a  single  arc 
(h,k)  .  Considering  all  nodes  (h)  with  c^  >  0  and  connecting  them  to 
all  k  such  k  e  D(h)  and  c^  <  0  we  have  constructed  a  bipartite 


66 


2 

graph  with  the  nodes  corresponding .  to  positive  as  one  set  of  nodes  and 
the  nodes  with  nonpositive  as  the  other  set.  Note  that  all  arcs  lead 

from  the  positive  set  to  the  nonpositive  set.  Figure  3.3  Illustrates  the 
original  network  and  its  bipartite  equivalent. 

The  structure  of  the  equivalent  bipartite  graph  of  the  network 
corresponding  to  the  dual  problem  (3.3.2)  resembles  that  of  a  simple  trans¬ 
portation  problem.  The  problem  could  be  solved  as  a  transportation  problem 
but  in  doing  so  by  the  primal-dual  method  of  Ford  and  Fulkerson  [9]  one 
essentially  only  need  solve  a  maximum  flow  problem  which  by  itself  seems  the 
most  attractive  method  of  solution. 

Consider  the  following  maximum  flow  problem  and  its  dual  where  we  assune 

*  °  *nd  i  0  • 


(3.3.3) 


Maximize  1*4 
i 


Subject  to 


0 

0 


Vi 

VJ 

V  i/c4  >  0 

V  J/ct  >  0 


2 

Bipartite  graph:  A  graph 
into  two  disjoint  sets  V  ■ 
(n^.nj)  c  A  joins  a  node 


G  *  (N,A)  in  which  the  node  set  N  decomposes 
jn^|  and  V  »  jn^j  such  that  each  edge 


n^  c  V  and  a  node  n^  e  V  . 


1 


(b)  Equivalent  Bipartite 

FIGURE  3.3:  EXAMPLE  OF  ORIGINAL  AND  EQUIVALENT  BIPARTITE  GRAPHS 
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The  dual  to  this  problem  is: 

Minimize  l  g.c.  -  l  he 

i  1  J  3  J 

Subject  to  -x^  +  >1  Vi 

(3.3.4)  x4  -  Xj  >0 

Xj  +  hj  .  0  V  j 

«i : 0  •  hj  >  o 

The  complimentary  slackness  condition;;  for  this  set  of  problems  are  for 
all  1  and  j  : 


fi<cl  “  fsi> 

-  0 

(a) 

hJ (_CJ  “  fjT) 

-  0 

(b) 

(§1  “  xi  *  1)f,i 

-  o 

(c) 

<X1  "  xJ)£lj 

-  0 

(d) 

(Xj  +  hj)fJT 

-  0 

(«) 

It  is  well  known  that  the  coefficient  matrix  of  problems  (3.3.3)  and 
(3.3.4)  is  unimodular  [2,  13]. 

The  unimodular  property  Insures  that  the  optimal  values  of  the  dual 
variables  are  Integers  but  not  that  they  are  only  H0"  or  +1  which  is 
required  to  relate  problem  (3.3.4)  and  (3.3.1).  The  following  theorem  will 
aid  In  establishing  the  zero-one  property  of  problem  (3.3.4). 

Theorem  3.3.1: 

If  .  (xj)  .  end  { hj  }  are  the  optimal  variable  a  for 

problem  (3.3.4),  then: 


» 


«  0  or  1  ;  x^  *  0  or  -X 


hj  ■  0  °r  1  I  *j  -  0  °r  -1 


Proof : 


For  feasibility  >  0  so  assume  g^  »  6  >  1  .  From  the  complimentary 
slackness  conditions  (3. 3.5a)  and  (3.3.5c),  this  implies  f^  ■  c^  >  0  and 
hence  xA  »  ggi  -  1  *  a  >  0  .  By  the  primal  feasibility  conditions,  we  must 
have  some  >  0  and  by  condition  (3.3.5d)  «=  x^  ■  a  >  0  .  Since 

f^  >  0  ,  it  must  be  that  fjT  >  0  which  makes  -h^  -  x^  *  a  >  0  by 
(3.3.5e)  so  that  h^  <  0  which  contradicts  the  feasibility  constraint 
hj  >  0  .  Thus,  0  <  g^  <  1  must  hold.  Since  -c^  *  0  has  no  effect  on  the 
problem  we  can  safely  assume  that  all  -c^  >  0  and  then  a  similar  argument 
goes  through  for  h^  showing  that,  0  <  hj  <  1  .  By  the  unlmodular  property, 
the  Kuhn-Hoffman  Lemma  given  in  Section  3.2.3  and  the  bounds  Just  shown  for 
gA  and  hj  we  have  that  in  the  optimal  solution  of  (3.3.4): 

-  0  or  1 


Also  we  have  that 


hj  -  0  or  1 


-x.  >  1  -  g  .  >  0  or  x.  <  0  , 
i  *  °si  *  j  ■  * 


xj :  ‘V  :  -1 


0  i  »i  ;  :  -i  • 
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Hence  -1  <  x.  <  0  and  -1  <  x.  <  0  . 

*  i  ■  ■  J  ■ 

x^  ■  0  or  -1  and  x ^  ■  0  or  -1  . 


Then  similarly  to  and 


I 


Q.E.O. 


Theorem  3.3.2: 

The  optimal  solution  to  the  subproblem  (3.3.1)  may  be  obtained  from  the 
optimal  solution  of  (3.3.4)  by  letting: 

\  • 1  -  «i 


Proof : 


Making  the  subsituation  as  given  above  for  g^  and  In  (3.3.4)  and 
using  the  bounds  established  In  Theorem  3.3.1  yields: 


Minimise 

Subject  to 
(3.3.6) 


}  \  -  l  Vi 

l/ci>0  i/c^O 


® :  \ : 1 
0  <  Yj  <  1 


Z  Vj  ■ « 


j/Cj<0 


>  0 
>  0 


+  V° 


(e) 

(b) 

(c) 

(d) 

(e) 

(f) 


Since  the  objective  functions  of  (3.3.6)  end  (3.3.1)  are  equivalent,  ve 

i 

need  only  establish  the  equivalent  feasibility  of  the  respective  optimum 

« 

solutions.  Because  (3.3.1c)  and  (3.3.6e,f)  are  the  same,  only  the 
equivalent  feasibility  of  (3.3.1b)  and  (3.3.6b,c,d)  need  be  considered. 

Let  D(h)  »  D~(h)  U  D+(h)  with  respect  to  the  dual  network  of  (3.3.1). 


Where: 


I 
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D  00  *  (j/h  >  j  and  <  0} 

and 

D+(h)  ■  { i/h  >  1  and  >  0}  . 

Every  feasible  solution,  hence  every  optimum  solution,  of  (3.3.1)  la 
feasible  for  (3.3.6)  since  (3.3.1b)  requires  that  if  Y^  »  1  then  •  1 
for  all  k  e  D(h)  and  all  that  (3.3.6b,c,d)  requires  Is  •  1  for  all 
j  e  D+(h)  CD(h)  . 

To  show  every  optimum  solution  of  (3.3.6)  is  feasible  for  3.3.1,  we 
observe  that  the  conditions  (3.3.6b,c,d)  imply  that: 

Yi  -  l  .>  *i  -  -l  »  -j  -  -1  ■>  r}  ■  1  v  i  c  d-(d 

This  leaves  only  to  establish  that  if  Y^  *  1  then  Y^  *  1  for 
k  e  D+(h)  .  If  this  condition  were  not  true,  then  the  value  of  the  objective 
function  (3.3.6a)  could  be  decreased,  since  c^  >  0  ,  contradicting  the 
optimality  assumption. 

Q.E.D. 


It  has  now  been  established  that  solving  the  maximum  flow  problem 

(3.3.3)  and  using  the  dual  variables  thus  generated  will  provide  a  solution 

to  the  subproblem  (3.3.1).  An  efficient  method  for  solving  problem  (3.3.3) 

is  essentially  the  labeling  procedure  developed  by  Ford  and  Fulkerson  [9,  21]. 

From  the  complimentary  slackness  conditions  (3.3.5),  for  problems  (3.3.3) 

and  (3.3.4),  we  know  that  if  it  is  impossible  to  saturate  any  arc  (s,l)  , 

that  is  when  f  ^  in  the  optimal  solution,  then  g^  *  0  and  hence 

Y^  ■  1  from  Theorem  3.3.2.  Also  if  some  f^  <  -c^  then  h^  ■  0  and 

Y.  *  0  .  From  Theorem  3.2,  we  know  that  if  Y.  *  1  ,  then  Y.  ■  1 
j  i  k 


for  all 


k  e  D(i)  .  It  also  follow*  that  if.  f^T  <  -c^  implying  that  ■  0  then 

Y.  ■  0  for  all  h  such  that  j  c  D(h)  .  These  conditions  are  key  factors 
n 

in  the  algorithm  to  be  discussed  and  presented  in  the  following  sections. 


I 
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3.4  PROPOSED  METHOD  OF  SOLUTION 

The  proposed  method  for  solving  the  original  subproblem  (3. 2. 1.2)  la  to 

transform  it  into  problem  (3.3.1);  transform  this  into  problem  (3.3.3)  aftd 

solve  (3.3.3)  as  a  maximum  flow  problem.  Since  the  network  corresponding  to 

(3.3.3)  is  bipartite,  actual  labeling  of  the  nodes  is  unnecessary  and  a  flow 

equal  to  min  (c^-c^)  for  k  >  j  can  be  allocated  at  each  Iteration.  If 

It  Is  impossible  to  route  all  flow  away  from  a  node  when  c^  >  0  then 

Y.  ■  1  and  from  Theorem  3.2  all  Y  ■  1  for  p  c  D(k)  .  Once  a  set 
t  P 

Y  -  (k/Y^  *  1)  with  fik  -  0  for  k  c  Y  ,  i  f  Y  has  been  determined  It  may 
be  removed  from  all  further  consideration.  The  values  of  must  be 

recorded  since  re-routing  of  flow  may  be  required.  Re-routing  may  be 
necessary  if  for  some  c^  the  following  conditions  hold: 

ck  >  l  m  ("Cj  "  [  flJ 
JcD(k)0S  1 

and  f  >  0  for  1  i  D(k)  ,  j  e  D(i)  fi  D(k)  where  S  ■  {j/c^  j*.  0)  • 

The  transformation  from  (3.3.1)  to  (3.3.3)  is  depicted  in  a  network 
sense  in  Figure  3.4.  Although  the  network  for  (3.3.3)  looks  like  a  maze, 
hence  making  the  problem  difficult  to  deal  with  in  this  form,  this  is  not  the 
case,  since  given  any  node  number  and  the  allowable  mining  sequence,  we 
know  exactly  which  nodes  to  look  at.  In  this  sense,  the  solution  procedure 
deals  with  the  network  in  Figure  3.4a  in  a  node-path  manner. 

The  dual  problem  (3.3.2)  variables  (p^)  are  just  p^  ■  0  for 
Cj  <_  0  and  ■  c^  -  f  otherwise.  In  effect,  we  have  solved  problem 
(3.3.2)  neglecting  the  {p^}  variables  and  determined  them  from  the  slack 
in  the  equations  corresponding  to  c^  >  0  . 
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(•)  Example  (3.3.1)  Network  (All  Unattached  Arcs  Leading  Out 
Go  to  Sink;  Pointing  into  a  Node  Are  From  Source) 


SOURC 


(b)  Example  (3.3.3)  Network  [The  (0,kij)  on  the  Arcs  Indicate 
the  Lower  and  Upper  Flow  Bounds] 

FIGURE  3.4:  EXAMPLE  NETWORK  FORMS  FOR  PROBLEMS  (3.3.1)  AND  (3.3.3) 
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The  other  dual  variables  {p^}  may  be  obtained  fiom  the  flow 
variables  {f^}  in  problem  (3.2.2)  as  follows: 


f^  *  flow  in  path  from  i  to  j  in  both  networks. 


Consider  any  arc  (p,q)  in  a  path  from  i  to  j  .  Set  p  ■  f 

if  (p,q)  is  not  in  any  previously  considered  path;  p  ■  p  +  f  If 

pq  pq  lj 

(p,q)  is  in  a  previously  considered  path. 

An  outline  of  the  subproblem  algorithm  developed  in  the  above 
discussion  follows. 


* 
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3.5  SUMMARY  OF  SUBPROBLEM  ALGORITHM 
Step  1.  Initialize: 

Let  S+  ■  {i/ct  >  0) 

S'  -  (J/Cj  <  0} 

Y  -  (Yk)  -  0 

Step  2.  Scan  Node  Min  k  e  S+  ! 

If  S+  ■  0  ,  go  to  Step  5. 

If  ck<  l  -Cj  •  V(k)  ,  go  to  Step  3. 
jeDOons’ 

Otherwise,  Step  4. 

Step  3.  Allocate  c^  : 

(a)  Let  ■  max  -  Cj 

JeD(k)ns” 

(b)  Set  ffcl  -  min  (c^.Cj) 

c^  ■  c^  -  and  delete  i  from  D(k) 

c.  ■  c,  -  {,. 
k  k  kl 

(c)  When  c^  ■  0  delete  (k)  from  S+  and  go  to  Step  2. 
Otherwise,  3(a). 

Step  4.  Determine  Profitable  Blocks: 

(a)  If  c.  >  V(k)  +  I  f..  -  V(k) 

*  i^D(k) 

JcD(i)pD(k) 

Set  Y  •  1  for  p  e  D(k)  and  delete  D(k)  from  S+  U  S 
P 

c^  -  c^  +  f y  for  i  {  D(k)  ,  J  {  D(k)  and  repeat 
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Steps  2-4  for  such  i  when  f ^  >  0  .  Then  go  to 
Step  2. 

(b)  If  ck  <  V(k)  ,  let  Q  -  {i/i  i  D(k)  ,  f  >  0  ,  j  c  D(i)  n  D(k)}  . 

For  each  i  e  Q  and  all  j  e  D(i)  f|  D(k)  with  f  >  0  ,  let 

4f^j  -  min  ,  AfRj  *■  min  (ck,-(Cj  -  Af ) )  , 

fkj '  fkj  +  s£kj  an,i  \  ■  'k  ■  A£kj  •  Thcn  vlth  £ij  ■  £u  - 
Af^  and  <=  Af^  repeat  Step  2.  Terminate  examination  of  Q 

if  ck  *  0  and  go  to  3(c).  If  Q  ■  0  and  ck  >  0  go  to 

Step  4(a). 


Step  5.  Terminate: 


Y  ■  { Y.  )  solves  the  subproblem  with  value  J  —  Y  *  £  c.  «  Z  , 

o  _  _  k  Ck  keS+  k 

where  —  *-  original  c.  ,  and  c.  *  modified  c.»  c.  -  T  f.  .  for 

<-k  K  k  K  k  j  kj 

k  £  S+  *  {k/^  >  0)  . 


It  is  quite  obvious  the  algorithm  terminates  in  a  finite  number  of 
steps  since  S+  is  finite  and  each  iteration  through  Steps  1-4  eliminates  at 
least  one  element  from  S+  and  the  procedure  is  stopped  when  S+  *  0  . 

Now  that  a  rather  simple  method  of  solving  the  subproblem  has  been 
developed  the  entire  problem  (3.1.1)  can  be  solved  by  the  decomposition 
principle  as  described  in  Section  3.2.1.  In  actual  computational  scheme,  it 
may  be  wise  to  terminate  the  subproblem  algorithm  once  Z  >  S  and  Introduce 
the  corresponding  { Y>  into  the  master  basis  even  though  the  subproblem  is 
not  optimal.  If  the  set  S+  is  large,  this  is  especially  recommended. 

A  general  flow  chart  of  the  subproblem  algorithm  follows. 
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GENERAL  FLOW  CHART  OF  SUB-PROBLEM  ALGORITHM 


Not  now! 


Remove  block  (k) 
from  set  of  posi¬ 
tive  value  blocks. 


Select  some  block  (l)€Q. 
for  each  block  (j)  in 
both  D(k)  and  0(1)  with 
fij  >  0,  allocate  profits 
from  (k)  to  (j)  and  from 
J  to  I  ,  then  to  some 
other  block  in  0(1)  if 
feasible. 


Have  alt  blocks 
with  positive  value 
been  scanned? _ 


Scan  remaining  positive 
block  (k),  nearest  surface: 
Is  it  profitable  to  mine 
block  (k)? 


Define  blocks,  (i)  ■  Q, 
not  in  D(k)  from  which 
profits  have  been  al¬ 
located  to  blocks  fj}  , 
in  D(k)  . 


Terminate, 
have  solution 


Classify  block  (k) 
and  its  restricting 
set  (D(k) )  as 
mineable.  Remove 
these  blocks  from 
further  consideration. 


If  profits  were 
allocated  to  D(k) 
from  other  blocks 
outside  this  set 
redistribute  such 
values  to  original 
blocks  and  restore 
original  blocks  to 
ositive  value  set. 


/Is  it  possible  to  al- 
Y es/  locate  all  profits  from 
\block  (k)  by  rerout inq? 
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3.6  NUMERICAL  EXAMPLE  OF  SUBPKOBLFM 

As  an  example  to  demonstrate  the  algoritlun,  consider  the  problem  In 
Figure  3.4. 


c,  *  4  <_  -c.  -  c_  *  5  .  Therefore,  go  to  Step  3. 


Therefore,  go  to  Step  4 


80 


Step  4a: 

■  5  >  3  +  0  .  The  last  term  is  aero  since  all  f^  ■  0  for  all 
J  c  D(ll)  and  i  i  D(ll)  .  Hence,  Yu  “  *  »  "  *2  "  Y3  "  Y6  "  Y7  "  1 

and  these  nodes  can  be  deleted  from  S+  and  s”  .  At  this  point,  the 
problem  has  been  reduced  to: 


s+ 

Ji 

s" 

12 

4 

4 

3 

13 

1 

5 

4  - 

14 

1 

8 

1 

15 

1 

9 

2 

Since  S+  la  not  empty,  return  to  Step  2  and  scan  node  12. 
Step  2: 

£  -CA  +  (-Cg)  ■  4  .  Therefore,  go  to  Step  3. 

Step  3: 

f12,4  "  3  *  f12,8  "  1  »  c12  -  0  ,  -cA  -  0  ,  -cg  -  0  . 

Step  2: 

c13  ■  1  <  -c4  +  (-Cj)  +  (-Cg)  +  (-Cg)  -  5  ;  go  to  Step  3. 
Step  3: 

j  ■  1  ,  -Cj  ■  3  ,  lu  ■  0  ;  others  unchanged. 

Step  2: 

•  1  1  -c4  +  (-c5)  +  (-Cg)  +  (-c9)  -  5  ;  go  to  Step  3. 
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Step  3: 


f14.5  '  1  >  -c5  ‘  2  •  c14  ’  0 


Step  2: 

c15  ■  1  <  -c4  +  (-c5)  +  (-cg)  *=  A  ;  go  to  Step  3. 


Step  3: 

f15,5  “  1  *  _C5  "  1  ’  Ci5  "  0  •  Now  all  nodes  of  S+  have  been 
examined  so  that  S*  is  considered  empty  so  go  to  Step  5. 

Step  5.  Terminate: 

The  following  {Y^}  give  the  optimal  solution  with  value  max  Z 


original 

c. 


final 


S" 

original 

“'I 

l 

Yi 

final 

- 12 

1 

3 

1 

0 

2 

2 

1 

0 

3 

1 

1 

0 

A 

3 

0 

0 

5 

5 

0 

1 

7 

1 

1 

0 

8 

1 

0 

0 

9 

2 

0 

2 

NOTE:  l  ljYi  +  [  CjYj  -  ^  ci  (final) 


max  Z  . 


The  values  of  the  flow  variables  are: 


f61  "  3  *  f62  "  1  *  f 10, 5  "  1  *  fll,2  "  1  *  f 11,3  "  1  »  f 11, 7  *  1  * 
f12,4  *  3  ’  f 12,8  ’  1  ’  f 13, 5  *  1  ’  flA,5  *  1  and  f15,5  "  1  * 


Hence,  the  are: 


I 


I 
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U61  •  3  •  “62  ’  2  •  u73  ■  1  1  U84  ■  3  •  “95  ■  2  •  °10,5  ’  2 
uU.6  ’  1  •  “11.7  ’  2  •  “l2.8  ‘  *  •  “13.9  '  1  ’  “U.9  '  1  1 
u35.10  '  1  • 

and  all  others  zero. 
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3.7  MODIFICATION  OF  PROCEDURE  WHEN  ASSUMPTION  2  IS  RELAXED 

Assumption  2,  as  discussed  in  Chapter  2,  constrained  the  blocks  of  the 
problem  to  be  all  of  equal  size.  The  result  of  this  was  that  the 
coefficients  of  the  dual  subproblem  variables  {p^}  (see  problem  3. 2. 2. 2) 
were  all  ones. 

The  method  of  solution  proposed  in  that  last  section  actually  neglected 
the  set  of  variables  { p ^ }  .  This  indicates  that  a  similar  procedure  will 
work  when  the  block  size  differs  in  certain  areas  of  the  deposit.  It  may  be 
desirable  to  have  different  block  sizes  in  various  zones  of  a  pit  due  to 
physical  elements  such  as  varying  pit  slopes. 

Consider,  for  example,  the  case  with  two  block  sizes  and  B2 

(volumes).  In  this  situation,  the  coefficients  of  the  variables  (p^ } 
would  be  Bj  and  B2  .  From  a  network  viewpoint,  we  then  have  two  sets  of 


nodes  defined  as: 


■  {j/block  (j)  has  volume  B^} 

N2  *  (j/block  (j)  has  volume  B2>  . 


The  maximum  flow  algorithm  of  Section  3. A  can  now  be  used  to  solve  the 
problem.  First  apply  the  algorithm  to  the  set 

U  (k/k  e  N2  ,  D(k)  p  t  0}  (assume  B^  <  Bj)  .  This  yields 

Y1  "  Y1  U  Y2  ‘  Where 


yJ  -  (y  /y  -  1  ,  D(j)  n  n2 


Y2  -  (Y  /Y  -  1  ,  D(j)  p  N2  /  0)  . 


and 
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CHAPTER  .4 

GENERAL  T-PERIOD  PROBLEM 


4.1  INTRODUCTION 

Open  pit  mine  production  scheduling  as  defined  in  Chapter  1  is  a 
very  real  problem  that  concerns  most  mine  operating  management.  The  best 
schedule  for  the  entire  planning  horizon  of  T-periods  is  usually  not  attain* 
able  by  repetitive  single-period  scheduling  because  of  the  dependence  of 
the  future  on  the  past  and  present.  Thus,  the  effect  of  any  single-period 
schedule  on  future  period  schedules  must  be  considered  in  a  multi-period 
scheduling  problem.  The  following  formulation,  developed  in  Chapter  2, 
illustrates  this  dependence. 

1  2  T 

Maximize  c^X  +  c2X  +  .  .  .  +  C.J.X 


(4.1.1) 


OjX1  +  02X2  +  .  . 

.  +  dtxt 

*  d 

AjX1 

“b1 

A2X2 

*  b2 

•  a/ 

bT 

ex' 

S  0 

BX1  +  BX2 

§  0 

1  2  3 

BX  +  BX  +  BX^  + 

.  .  .  +  bxt 

§  0 

£  Xql  ♦  E  xq2  + 
h  h 

4  q 

.  .  +  £  XqT 
h 

q 

s  1  ¥  h 

xq  s  0 

*h  "  u 
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Host  of  the  concepts  in  this  chapter  will  be  examined  in  the  context 
of  a  two-period  problem.  General  extent  ions  will  be  presented  only  when 


they  do  not 

follow  direct ly. 

Thus  the  probl 

em  to  be 

considered  has 

form; 

Maximize 

CjX1  +  c2X2 

(•> 

Subject  to 

DjX1  +  02X2 

■  d 

(b) 

AjX1 

*  b1 

(c) 

a2x2 

"  b2 

(<0 

(4.1.2) 

BX1 

*  0 

(e) 

BX1  +  BX2 

&  0 

<0 

E  xjj1  +  E  Xj2 
q  q 

*  1 

(9) 

X?*  *  0 
h 

(h) 

In  any  practical  situation  this  problem,  as  in  the  case  of  the  prob¬ 
lem  discussed  in  Chapter  3,  rapidly  exceeds  the  capabilities  of  available 
linear  programming  codes,  even  with  the  most  modern  computing  facilities, 
and  a  more  compact  method  of  solution  is  necessary.  Because  of  the  special 
structure  of  the  dynamic  constraints  (4.1.2),  and  especially  of  the  B 
matrices,  some  form  of  decomposition  seems  promising. 

4.2  DECOMPOSITION  OF  GENERAL  PROBLEM 

There  are  many  ways  to  decompose  (4.1,2).  The  problem  as  formulated 
has  essentially  two  parts:  (1)  the  blending  and  allocation  portion  corres¬ 
ponding  to  the  { Dh}  and  |A^|  ,  and  (2)  the  allowable  mining  sequence 
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portion  characterized  by  the  B  matrices.  The  most  attractive  decomposi¬ 
tion  approaches  to  the  problem  appear  to  be  those  which  permit  the  reduc¬ 
tion  of  the  B-structure  by  virtue  of  Theorem  2.1,  and  allow  isolation  of  • 
the  EY*1  5  0  as  sub-problems,  hence  using  their  simplified  structure  to 
advantage.  It  is  reasonable  to  expect  that  the  blending  and  allocation 
part  of  the  problem  [4.1.2  (b)  -  (d)]  ,  will  be  of  such  a  dimension  that 
It  can  easily  be  solved  as  a  linear  program,  either  as  a  master  or  sub-master. 

This  rationale  leads  to  the  following  master  problem: 
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so  that  the  sub-problem  becomes; 


(if. 2. 2) 


Maximize  c^X*  +  c^X^ 

Subject  to  BX1  §  0 

BX1  +  BX2  *0 

E  x'1'  +  £  X^2  s  )  V  h 


?  0 


By  virtue  of  Theorem  2.1  this  can  be  reduced  to 


Maximize  c^Y*  +  (a) 

Subject  to  EY1  SO  (b) 

EY1  +  EY2  S  0  (c) 

(4. 2. 3)  IY1  +  IY2  S  1  (d) 

Y*  £  0  (e) 


Where  ct  *  {^s t }  suc^  that  St  "  max  St  an^  Y*  *  {Yh  }  such  that 
Yh  *  Xh^  wMch  corresponds  to  6^  . 

The  problem  as  given  in  (if. 2.3)  has  many  less  variables  than  (4.2.2) 
and  has  a  structure  that  is  much  easier  to  deal  with.  This  reduction  step 
amounts  to  a  classification  and  selection  of  the  method  of  treatment  of  a 
block  If  it  is  to  be  mined.  It  reduces  thz  sub-problem  strictly  to  a 
scheduling  problem  in  which  the  remaining  decisions  are  if  and  when  to 
select  a  block  for  mining. 


* 
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4.3  SUB-PROBLEM  SOLUTION 

The  sub-problem  (4.2.3).  is  still  a  rather  large  problem  and  not  easily 
solved  in  a  straightforward  manner.  The  triangular  structure  could  be  cm* 
ployed  to  some  advantage  but  is  somewhat  cumbersome  due  to  the  dimension 
of  the  problem. 

By  making  the  proper  substitution  of  variables  (4.2.3)  can  be  trans- 

formed  Into  a  more  workable  form.  Let  w  *  E  Y  .  Then  (4.2.3) 

k  -  I 

becomes: 


Maximize 

(2,  -  «2)w' 

.  2 
+  c?w 

(a) 

Subject  to 

Ew1 

£  0 

(b) 

Ew2 

£  0 

(c) 

Iw2 

£  I 

.<“> 

,1  .2 

Iw  •  Iw 

£  0 

(®) 

t 

w 

£  0 

(0 

In  (4.3.1),  (e)  and  (f)  replace  (4.2.3)(e).  Thus  the  transformation 
has  increased  the  constraint  set,  but  the  problem  now  has  a  more  computa¬ 
tional  efficient  structure. 

In  Chapter  3  an  efficient  method  was  proposed  to  solve  a  problem  of 
the  form 

Maximize  cY 

Subject  to  EY  5  0 

(4.3.2)  IY  5  1 


Y  £  0 


I 
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which  used  the  structure  of  the  matrix  E  to  advantage.  It  seems  attrac¬ 
tive  to  again  isolate  problems  of  the  form  (4.3.2)  and  thus  use  the  method 
de  veloped  in  Chapter  3.  It  is  clear  from  the  structure  of  problem  (4.3.1) 
that  the  decomposition  principle  could  be  applied  again  in  a  straightforward 
manner  with  the  constraints  (4.3. 1  e )  as  the  master  problem  with  P  sub¬ 
problems  each  with  a  form  corresponding  to  (4.3.2).  In  such  a  scheme  the 
master  problem  wou Id  be  relatively  large  and  even  though  the  constraints 
are  simple  it  would  require  sufficient  calculations  to  introduce  the  sub¬ 
problem  solutions  into  the  master  basis  to  make  this  suggestion  unattractive. 

It  is  known  from  the  Theorem  of  Hel ler-Thompkins-Ga le  given  in  [13] 
that  the  coefficient  matrix  for  problem  (4.3.1)  is  unimodular  and  hence  from 
the  Hoffman-Kuhn  Lemma,  [13],  stated  in  section  3.2.3  that  the  optimal 
{wj}  are  zeros  and  ones.  This  fact  and  the  implications  of  the  constraint 
set  (4.3.  le)  lead  to  an  efficient  method  of  solution  through  decomposing 
by  both  rows  and  columns.  The  constraint  set  (4.3.  le)  implies  that  problem 

(4.3.1)  may  be  solved  by  considering  a  sequence  of  problems  in  the  form  of 

k  t  h 

(4.3.2)  but  need  only  consider  those  w^  in  the  k  stage  problem  which 

k  +  1 

correspond  to  w^  *  I.  To  solve  problem  (4.3.1),  first  consider; 

2 

Maximize  CjW 


2 

Subject  to  Ew  £  0 

(4.3.3)  Iw2  5  0 

w2  5  0 


Solving  this  problem  yields  an  optimal  solution  w2 


2  2 
W*  U  VT 
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where 

and 

Next,  consider 


(*♦.3.4) 

where 

.1  )  1  .  -2 

w  .  |wh  /  wh 

The  optimal  solution  of  this  yields  w  *  w|  U  W*  ,  with  w|cyj  . 

In  the  terms  of  a  mining  schedule,  problem  4.3.3)  determines  those 
blocks  which  are  profitable  for  mining  in  periods  I  and  2  ;  then  solving 
problem  (4.3.4)  selects  from  this  group  those  which  are  best  taken  in 
period  I  . 

Next  we  must  consider  blocks  which  were  not  profitable  for  mining  in 
period  2  but  whose  value  in  period  I  may  be  such  that  they  can  be  pro* 
fitably  taken  in  period  I  .  This  may  also  force  blocks,  which  have  already 
been  selected  for  period  2  ,  to  be  rescheduled  into  period  I  . 

To  determine  such  blocks  consider  only  the  set  W*  and  adjust  the  cost 


vl  ■  {"h  '  ;,h  - '} 

w?"{;h  '"i;  *°j 


Maximize  (c^  -  c^Jw' 


ca1 

Ew 


Iw 

.1 

w 


I 


£  0 


5  I 


£  0 


- 


coefficients  of  these  variables  as  follows: 


if  h  \  w  n  w' 


if  h  t  w1  n  w 2 


Solve 


Maximize  c.w 


Subject  to  Ew  r5  0 


Iw  =E  i 


w  *  0 


- .  112 
giving  the  optimal  solution  w‘  ,  Then  alter  the  sets  W  ,  W  ,  W 


and  W  as  fol lows; 


w!  *  jwh '  \  *  '•  “h 

w- '  {wh 7  w‘h'  *  °| 

vl  ‘(wh '  “h  ■  '• 1 


h  ■  '•  wh 


J.  t  2  .  -I 

W-  ‘  K  1  wh 


These  sets  provide  the  optimal  solution  to  problem  (4.3.1). 

The  steps  of  the  algorithm  will  now  be  presented  in  terms  of  the 
general  T-period  sub-problem  (4.3.5)  where  c  *  c,  -  c 

K  K  K  +  I 
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Maximize 

-  2 
C1W 

-  2 

+  C,jW  +  .  . 

.  + 

Subject  to 

Ew1 

S  0 

• 

N 

5 

Ui 

S  0 

• 

*  EwT 

S  0 

IwT 

S  1 

Iw1  - 

Iw2 

S  0 

Iw2  -  Iw3 

S  0 

•IwT"1  -  IwT  SO 


wl  ?  0 
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OUTLINE  or  SU6- PROM.EH  ALGORITHM 

STEP  I.  INITIALIZE: 

Set ,  Vk  =  |a  1 1  w^} 

T  =  T 

Yk  =  ()  ,  Yk  =  |all  wk  |  for  all  k  . 

STEP  2.  SOLVE  STAGE  k  PROBLEM: 

,  .  .  -  k 

Maximize  c.  w 

k  k 

V 

le 

Subject  to  Ew  SO 
Iwk  S  1 
wk  £  0 

Let  w  be  the  optimal  solution  and  define; 


Add  Wk  to  Yk  and  delete  from  Yk  and  if  k  =  T  add  Wk  to 
Y*  and  delete  from  Yl  for  t  *  k  +  1 . T. 

STEP  3. 

If  T  *  I,  terminate.  The  sets  Yk  and  Yk  for  k  *  1,  2,  ...  T 
solve  the  sub-problem,  if  k  =  I,  set  k  =  T  and  go  to  STEP  4. 
If  Vk  '  #  ()  ,  set  k  =  k  -  1  and  go  to  STEP  2.  Otherwise  go 


to  STEP  4. 


I 
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STEP  4.  ADJUST  Cf  _  , 

Set 

.  I  "  -  |  +  for  h  €  Yf 

Set 

k  -  k  -  I 


T  ■  T  -  1  and  go  to  STEP  2. 


THEOREM  4.1:  T/ie  tub-ptioblm  algorithm  temninalu  in  a  friniXe.  numbin 
oh  4tip*  uiith  thi  optimal  solution  oh  problem  (4.3.5). 

Proof:  Since  the  algorithm  starts  with  period  T  (T  ■  T,  k  ■  T)  and  pro¬ 
ceeds  step  by  step  through  period  i  (T  «  1)  it  indeed  terminates  in  a 
finite  number  of  steps. 

Because  the  conditions  of  the  algorithm  always  insure  feasibility  of 
problem  (4.3.5) .  to  demonstrate  optimality  it  must  be  shown  that  there  does 
not  exist  a  node  (h)  such  that 

(a)  wk  ■  0  ,  k  ■  maximum  t  /  w*  ■  0  and  c|*  ♦  E  c!  >  0 

h  h  k  I  6  0(h)  k 

(b)  w*  «  1,  z  *  minimum  t  /  w*  ■  I  and  h  €  D(i)  /  c*  ♦  E  C^<0 

h  h  '  q  €  0(1) 


or 


% 


I 


Clearly  the  algorithm  does  not  allow  all  these  conditions  to  hold.  For 

k 

example,  consider  the  case  where  T  c  2  .  When  w,  =  0  and  k  e  2  then 


the  initial  pass  through  STEP  2  prohibits  condition  (a)  by  insuring  w.  *  I 


If 


i  €  0(h) 
£ 


'2  >  0  •  If  "h 


*  0  and  k  =  1  ,  then  by  the  algorithm 


lht  +  £  c!  <  0  ,  thus  condition  (a)  cannot  hold.  When  w*  *  I  and 

1  I  €  0(h)  1  b 


z  «=  2  ,  then  as  before  the  initial  STEP  2  assures  h  C  D(i)  such  that 
l!  +  £  c?  >  0.  If  wf  >=  I  and  z  «  1  ,  then  either  both 

n  /-  r./  2  h 

q  €  0(i) 

£*  ♦  £  cq  >  0  for  t  *  1  and  2  or  for  (i)  €  Y,  such  that 

1  q  €  D(  i)  1 

h  €  0(1)  ,  ci  +  £  cjsO  and  (i)  for  period  1  is  such  that 

q  €  0(i)  U  Y2 

cj  ♦  2  cq  -  cj  >  0  so  that  cj  -  *  cj  >  0. 

'  q  €  0(1)  11  '  2  q  €  D(i)  U  ^  Z  1 

Hence  condition  (b)  cannot  hold  and  Theorem  is  proven.  QED 


4.4  NUMERICAL  EXAMPLE  OF  2-PERI0D  SUB-PROBLEM 

Consider  the  following  networks  corresponding  to  a  2-period  sub-problem 
in  the  form  of  (4.2.3).  The  number  in  the  lower  half  of  the  circle  repre¬ 
sents  the  |c^  |  while  the  node  number  is  shown  in  the  upper  half. 


(a)  Period  1  (b)  Period  2 


% 


FIGURE  4.4.1; 


EXAMPLE  NETWORKS  FOR  TWO-PERIOD  PROBLEM 
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FIGURE  4.4.3:  NETWORK  FOR  PROBLEM  IN  FIGURE  4.4.2 
WHEN  T  *=  2  AND  k  »  1 


On  this  pass  STEP  2  yields 


and 


Since  k  *  i  in  STEP  3  proceed  to  STEP  4  which  yields  the 
in  Figure  4.4.4  with  T  *  I  and  k  =  1  . 


FIGURE  4.4.4;  EXAMPLE  NETWORK  FOR  T  »  I  AND  k 
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Another  STEP 


and  since  T 


In  other 
blocks  1-7,  9 


2  provides 


W 


.  1  1 
t w  I »  ^  * 


1  1  I  1  |  | 

V  V  V  w9’  wloJ 


"  *  »  termination  is  reached  in  STEP  3  with 

V  ■  |"|.  wj,  w^,  „j,  wj0| 

|we,  wM,  wj2j 
{al.  4} 


words  the  optimum  mining  plan  for  this  sub-problem  Is  to  take 
10  in  period  1  and  blocks  8,  II,  12  In  period  2. 


I 
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4.5  GENERAL  COMMENTS  AND  PROCEDURAL  SUMMARY 

The  type  of  decomposition  utilized  in  the  algorithm  presented  in  sec¬ 
tion  4.3  reduces  the  size  of  the  problems  with  the  form  (4.3.2)  which  need 
to  be  solved  at  any  one  time  by  the  methods  of  Chapter  3.  In  this  sense 
it  reduces  the  average  number  of  paths  incident  out  from  any  positive  node 
and  hence  should  decrease  computation  time  considerable.  Even  though  com¬ 
putational  experience  is  not  available  at  present  to  verify  this  assertion, 
there  seems  to  be  a  strong  indication  that  this  is  true.  In  [10],  Gilbert 
deals  with  a  comparable  scheme  and  his  conclusion  is  (from  the  viewpoint 
of  this  paper)  that  the  ability  to  partition  the  problem  (4.3.2)  into 
separate  parts  would  be  computationally  efficient.  Even  though  he  is  not 
actually  dealing  with  the  same  form  of  the  problem,  his  conclusions  seem 
applicable  here. 

The  extension  of  the  methods  described  in  section  4.3  to  the  case  of 
unequal  block  size  can  be  accomplished  with  ease  in  a  manner  similar  to 
that  discussed  in  Chapter  3. 

With  the  developments  of  Chapters  3  and  4,  the  T-period  open  pit  mine 
production  scheduling  problem  can  now  be  handles  quite  nicely  by  employing 
the  decomposition  principle  of  Datzig  and  Wolfe,  thus  taking  full  advantage 
of  efficient  network  techniques  permitted  by  the  simple  structure  of  the 
sub-problems. 

The  procedures  and  the  algorithms  developed  in  this  present  investi¬ 
gation  may  be  used  for  all  phases  of  mine  planning;  optimum  long  range  planning 
and  optimum  short  range  planning  as  well  as  optimum  operational  or  actual 
production  scheduling. 

In  order  to  utilize  the  methods  of  this  study,  one  must  first  of  all 
formulate  the  scheduling  problem  as  a  large  scale  linear  programming  problem 


considering  the  constraints  which  govern  the  system  as  discussed  In  Chapter 
2.  Next,  the  decomposition  principle  is  used  to  decompose  the  problem  Into 
a  rather  simple  linear  programming  problem,  called  the  master  problem,  and 
a  sub-problem  which  can  be  easily  solved  by  the  methods  developed  In  Chapters 
3  and  4.  The  master  problem  is  first  solved  to  obtain  the  pA/tcea  t  n  and 
s  ,  from  which  the  profit  coefficients  for  the  sub-problem  are  determined. 
After  the  subproblem  is  solved  for  a  mining  plan  based  on  the  profit  coef¬ 
ficients,  the  value  of  this  plan,  cY  ,  (see  problem  4.2.3)  is  compared  to 
the  price,  s  ,  for  opt imab i 1  i ty.  If  cY  >  s  it  is  profitable  to  include 
the  present  plan  in  the  master  problem.  This  procedure  is  repeated  until 
the  optimality  condition,  cY  5  s  ,  Is  realized.  The  master  program  selects 
a  convex  combination  of  the  sub-problem  plans  which  satisfies  the  given  con¬ 
straints  and  maximizes  total  profits.  To  initiate  the  master  problem  an 
initial  feasible  mining  plan  is  required.  Such  a  plan  is  readily  available 
since  any  mining  scheme  which  satisfies  the  allowable  mining  sequence  con¬ 
straints  wi 1 1  do. 

To  aid  the  presentation  of  the  complete  procedure  the  following  flow 
chart  is  given. 

In  the  next  chapter  the  flexability  and  general  utilization  of  the 
proposed  methods  will  be  discussed  and  illustrated  through  numerical 
examples. 
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GENERAL  FLOW  CHART  OF  COMPLETE  PROCEDURE 


Yes 


£ 


Is  it  profitable  to 
consider  new  sub-problem 
plan  in  master?  (cV  >  s) 


Solve  stage  k 
problem  for  mineable 
and  non-mineable 
blocks  using  algorithm 
given  in  Chapter  3. 


Have  a  1 1 
/  periods  beer 
\  examined? 

»  1) 


Yes 


Have  al 1  blocks 
found  profitable  for 
period  T  been  scheduled 
;  »  1  or  '  1  m  P)? 
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CHAPTER  5 

MINE  PLANNING  AND  SYSTEM  MANAGEMENT 

5.1  INTRODUCTION 

Mine  planning  and  system  management  contribute  greatly  to  converting  a 
mineral  deposit  into  an  economical 1y  feasible  mining  operation.  If  manage- 
ment's  goal  of  maximum  profit  is  ever  to  be  attained,  scientific  decision 
tools  similar  to  those  developed  in  this  investigation  must  be  substituted 
for  present  mine  planning  practices. 

Besides  providing  that  open  pit  mining  schedule  which  maximizes  profit, 
the  model  and  the  method  of  solution  discussed  in  previous  chapters  yields 
other  benefits  as  well.  The  techniques  of  sensitivity  analysis  by  parametric 
variation  of  costs  or  constraints,  provides  much  valuable  information  to 
aid  management  In  evaluating  their  policies,  mining  techniques,  system  de¬ 
sign,  and  technological  changes.  In  a  sense,  solving  the  model  is  a  "lab 
experiment,"  in  which  many  different  possible  plans  resulting  from  various 
policies,  practices  and  technologies  may  be  rapidly  expanded  to  evaluate  the 
ensuing  patterns  of  development.  The  procedures  presented  in  this  study  may 
also  affect  some  of  the  traditional  concepts  used  in  mine  planning  such  as 
cut-off  grade  or  economic  cut-off  in  the  sense  they  are  presently  known  to 
be  applied  as  discussed  in  Chapter  1.  Also  one  can  view  the  entire  planning 
horizon  through  the  model  and  thus  consider  the  interaction  of  demands  of 
period  schedules.  This  insures  the  attainment  of  the  maximum  total  profit 
goal  if  a  feasible  plan  exists.  This  goal  of  maximum  total  profit  may  not 
be  attainable  by  sequential  sub-optimization:  maximizing  profits  for  each 
period  in  succession  subject  to  only  the  constraints  of  that  period. 

In  order  to  illustrate  some  of  these  concepts  some  simple  numerical 


examples  will  be  presented. 
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5.2  NUMERICAL  EXAMPLES 

5.2.1  Model  Setting 

A  two  dimensional  model  based  on  the  data  in  Table  5.1  will  serve  to 
illustrate  how  mining  programs  (schedules)  may  vary  with  management  objec¬ 
tives  and  the  constraints  induced  by  the  system  or  management  policies. 

The  system  considered  here  is  a  mine  and  a  concentrator.  It  is  as¬ 
sumed  that  the  concentrating  process  is  not  varied  to  change  output  but  that 
it  is  operating  within  its*  most  efficient  range.  The  valuable  constituent 
within  the  ore  will  not  be  designated,  but  if  the  reader  desires  he  may  think 
of  the  ore  as  an  iron-bearing  material.  The  assay  va 1 ues  are  then  per-cent 
iron. 

5.2.2  Example  1 

The  ultimate  pit  limit  obtained  by  using  the  method  of  Chapter  3  (with¬ 
out  constraints)  is  shown  in  Figure  5.1.  The  total  tninabte.  as  in¬ 

dicated  by  the  ultimate  contour  has  a  value  of  $232,320  (96.8  x  2400) 

(profit  before  taxes).  This  •’i&icAve  includes  9600  cubic  yards  waste  and 
110,400  tons  crude  ore  yielding  42,219  tons  concentrate  with  an  average 
grade  of  66.05. 

Note  that  within  the  ultimate  limit  there  are  blocks,  for  example  (1,4) 
and  (0,12),  which  are  designated  as  ore  even  when  their  profit  values  are 
negative.  These  blocks  are  such  that  their  value  as  an  ore  is  greater  than 
their  value  as  waste.  Such  ore  blocks  with  negative  values  will  only  be 
mined  if  the  underlying  blocks  cover  the  loss  of  the  negative  ore.  For  ex¬ 
ample  if  block  (1,4)  had  been  designated  waste,  it  would  have  been  uneconomical 
to  recover  block  (2,5),  and  hence  total  profits  would  have  been  reduced  to 
96.0  instead  of  96.8  .  (These  values  are  the  profit  index  values.  The  true 


Ultimate  limit  (doubt 
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profit  may  be  obtained  by  multiplying  by  2^00.)  In  all  further  discussion 
the  index  profit  values  will  be  used. 

Now  let  us  assume  that  we  are  going  to  mine  this  deposit  over  a  time 
horizon  of  three  periods.  Also  assume  that  the  restrictions  on  the  system 
are  only  on  mining  capacity  per  period.  These  restrictions  are  that  not 
more  than  8  block  units  may  be  taken  in  period  I  and  not  more  than  10 
in  each  of  the  remaining  periods. 

The  problem  now  is  to  plan  the  operation  such  that  the  restrictions 
are  not  violated,  but  so  that  we  maximize  total  profit.  It  is  known  from 
the  unrestricted  ultimate  plan  that  the  maximum  profit  realizable  from  this 
deposit  is  96.8.  If  this  can  be  obtained  by  a  mining  program  that  satisfies 
the  capacity  restrictions,  such  a  program  is  optimal. 

From  Figure  5.2  it  is  seen  that  there  are  at  least  three  different 
schedules  which  meet  these  criteria.  Any  of  these  plans  is  optimal  with 
respect  to  our  assumed  model  and  could  be  obtained  by  the  methods  proposed 
in  this  study.  The  optimal  solution  in  this  case  is  not  unique  and  a  post- 
optimal  analysis  of  the  model  would  provide  all  the  possible  variations  of 
the  optimal  solution. 

In  plan  1,  Figure  5.2a,  the  profit  values  for  periods  1  and  2  are  more 
nearly  equal  then  in  the  other  plans.  In  plan  2  the  concentrate  tonnages  for 
periods  1  and  2  are  approximately  the  same.  The  cumulative  schedules  for 
periods  1  and  2  are  identical  for  plan  2  and  the  plan  shown  in  Figure  5.2c, 
yet  for  the  individual  periods  the  mining  plans  differ. 

Figure  5.2c  shows  the  mining  schedule  which  maximizes  the  profit  in 
period  1  subject  to  the  mining  capacity  constraints  and  then  from  what  remains, 
provides  the  maximum  profit  in  period  2,  again  subject  to  the  capacity  con¬ 
straints.  The  remaining  material  within  the  ultimate  contour  is  recovered 
in  period  3  and  also  satisfies  the  capacity  constraints.  This  last  plan 
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illustrates  sequential  sub-opt  imizot ion  which  in  this  case  also  maximizes 

i 

total  profit.  This  situation  is  not  always  true  as  will  be  Indicated  in 
•  the  next  section. 

The  mining  programs  shown  in  Figure  5.2,  thus  illustrate  that  there  may 
be  a  number  of  feasible  plans  which  satisfy  the  overall  objective  of  maxi* 
mum  profit.  However,  as  also  is  shown,  there  may  be  other  objectives  or  con¬ 
siderations  such  as  to  equalize  production  among  periods  or  even  to  maximize 
profits  in  a  particular  period.  Many  times  these  considerations  may  be  due 
to  what  is  termed  the  intuitive.  mining  appAoac.li  or  the  a**  <>i  mining  . 

With  the  techniques  developed  in  this  investigation  many  of  these  practices 
may  be  studied  and  evaluated. 

5.2.3  Example  2 

In  order  to  illustrate,  to  a  greater  extent,  the  effect  and  cost  of 
system  constraints  and  management  policy  let  us  now  Impose  further  restric¬ 
tions  on  the  example  model.  Assume  together  with  the  mining  capacity  con¬ 
straints  of  the  previous  example  we  have  the  followii.j  restrictions  on  the 
concentrator  and  concentrator  product: 

1.  Available  concentrator  hours  for  period  1  2  240. 

2.  Available  concentrator  hours  for  periods  2  and  3  5  480. 

3.  Average  grade  of  concentrate  be  between  64.0  and  66.0  in  each 
period. 

The  objectives  will  be  to  maximize  total  profit. 

As  was  indicated  in  the  discussion  of  the  ultimate  contour,  the  average 
grade  of  the  total  ultimate  reserve  was  66.05,  hence  the  last  restriction 
may  be  somewhat  unrealistic.  But  for  the  sake  of  illustration  we  will  as¬ 
sume  this  is  the  result  of  some  management  policy.  Because  of  this  grade 


restrict  i  oo  .is  wt,  i  I  a  3  the  others  v.e  know  that  our  total  profit  may  be  less 
th.in  96. P  and  moo  of  the  following  or  both  arc  going  to  occur.  First  in 
order  to  meet  the  grade  restriction  of  66.0  it  may  be  necessary  to  blend  in 
some  waste  material  with  the  ore.  Another  way  to  possibly  achieve  this  limit 
is  to  extend  the  mining  operation  beyond  the  optimum  ultimate  contour.  Using 
the  methods  developed  in  the  present  investigation  the  best  way  to  meet  the 
given  restrictions  will  be  selected  from  these  alternatives. 

Shown  in  Figure  5.3  are  two  mining  plans  for  our  three  period  problem 
which  satisfy  the  constraints  on  mining  capacity,  concentrator  capacity  and 
management’s  grade  restrictions.  These  plans  were  obtained  by  using  the 
procedures  developed  in  Chapters  3  and  4.  To  start,  one  of  the  plans  given 
in  Figure  5.?,  was  used  as  an  initial  feasible  subproblem  solution,  which 
even  though  infeasible  for  the  master  problem  provided  initial  prices.  After 
a  few  iterations  of  the  master  problem  and  the  subproblem  routine,  plan  1  was 
obta i ned. 

In  plan  1  the  period  schedules  remain  within  the  ultimate  limit  but  some 
of  the  material,  blocks  (0,6)  and  (0,11),  which  was  previously  designated 
as  waste  has  been  utilized  as  ore.  The  total  profit  at  this  stage  of  the 
master  problem  is  96.715. 

Plan  2,  Figure  5.36,  is  the  result  of  another  iteration  through  the  sub¬ 
problem  routine  and  the  master  problem.  The  pricing  mechanism  of  the  master 
problem  has  indicated  that  profits  may  be  increased  to  96.756  by  remaining 
with  the  waste  designation  of  blocks  (0,6)  and  (0,11),  except  for  a  small 
amount  of  block  (0,6),  but  extending  the  mining  operation  beyond  the  optimum 
ul  t imate  1 imi t . 

Figure  5.4  illustrates  the  optimum  mining  schedule  for  the  three  period 
problem  which  was  obtained  after  one  more  iteration.  The  average  grade  in 
each  period  is  66.0.  This  plan  is  similar  to  plan  2,  but  note  now  that  all 


FIGURE  5.3:  FEASIBLE  NON-OPTIMAL  MINING  PLANS  FOR  RESTRICTED  MODEL 
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of  both  blocks  (0,6)  and  (0,11)  are  designated  as  waste  and  the  extension 
beyond  the  ultimate  limit  is  slightly  less  than  in  plan  2.  The  total  maximum 
profit  is  96.781  and  thus  some  restriction  has  caused  a  decrease  In  profit 
from  that  of  the  ultimate  limit.  Since  as  is  shown  in  Table  5.2  there  still 
Is  concentrator  (plant)  capacity  remaining  for  each  period,  the  conclusion 
Is  that  the  decrease  in  profit  Is  mainly  due  to  the  grade  restriction. 

The  previous  example  has  demonstrated  how  the  cost  and  effect  of  system 
constraints  and  management  policies  may  be  evaluated.  Even  though  the  given 
grade  restriction  may  not  have  been  realistic  considering  the  average  grade 
of  the  entire  deposit,  it  still  brings  out  the  value  of  such  an  analysis. 

Many  of  the  system  constraints  and  policies  are  unavoidable,  but  knowing  their 
cost  and  effect  may  be  a  tremendous  aid  to  management. 

Figure  5.5  shows  the  mining  plan  for  period  1  which  yields  the  maximum 
profit  subject  to  the  mining  capacity,  concentrator  capacity  and  grade  con¬ 
straints  of  the  previous  problem.  This  is  the  initial  phase  of  determining 
a  mining  plan  by  sequential  sub-optimization  in  that  the  constraints  for  the 
other  periods  were  ignored.  This  plan  for  period  I  differs  considerably  from 
that  of  the  period  I  plan  within  the  maximum  total  profit  plan.  (Compare 
Figures  5.4  and  5.5)  Note  that  in  the  plan  of  Figure  5.5  it  is  more  pro¬ 
fitable  totreatsome  waste  blocks  (0,6)  and  (0,11)  as  ore  rather  then  ex¬ 
end  the  mining  beyond  the  ultimate  contour. 

The  period  I  profit  in  the  sequential  sub-optimization  plan  (Figure  5*5) 
is  41.427  which  is  nearly  triple  that  of  the  period  I  plan  given  in  Figure 
5.4.  This  looks  quite  good  until  we  consider  the  mining  system  over  the 
entire  planning  horizon  of  three  periods  and  the  objective  of  total  maxi¬ 
mum  profit. 

First  of  all  we  must  recall  that  any  deviation  from  the  optimum  ultimate 
contour  plan  costs  money.  Since  the  average  grade  of  the  remaining  material 
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0,7 

.7936 

.  2064 

1.000 

0,8 

.8320 

.  1680 

.3990 

0,9 

.6224 

.3776 

1.0000 

0,  10 

1.0000 

1.0000 

0,11 

1.000 

— 

1.0000 

0,12 

.6523 

.3478 

0,13 

.6733 

.3267 

1,3 

— 

.1875 

1,4 

.0919 

.9081 

1,5 

1.000 

1,6 

1 .000 

1,7 

1.000 

1.8 

.2348 

.7652 

1.9 

.2348 

.7652 

1,10 

.0750 

1.000 

1.0000 

1,11 

.3267 

.6733 

1,12 

.0919 

.9081 

2,4 

.1875 

2,5 

1.0000 

2,6 

1.0000 

2.7 

.6285 

.3715 

2,8 

.6507 

.3493 

2,9 

.3267 

.6733 

2,10 

1.000 

2,H 

- - 

1.000 

Tota  l 

6.7304 

9.9430 

up 

6.650 

Pit  Hrs 

225.17 

403.92 

m 

240.0 

!. 1478  ore 
.8522  waste 


in  the  sequential  plan  is  65.94  we  would  expect  that  the  grade  restrictions 
could  be  met.  We  would  also  expect  to  stay  within  the  concentrator  capacity 
limitations.  However,  consider  the  mining  capacity  constraints.  As  indicated 
in  Figure  5.5,  a  total  of  6.65  blocks  were  taken  in  period  1.  This  leaves 
20.35  blocks  to  be  mined  in  the  remaining  two  periods  if  we  are  to  adhere  to 
the  optimum  ultimate  plan.  The  condition  could  be  relieved  by  taking  some 
waste  blocks  in  period  1.  This  would  of  course  decrease  period  1  profit. 

No  further  ore  material  could  be  treated  in  period  1  since  the  concentrator 
is  at  maximum  capacity.  By  shifting  waste  blocks  we  would  not  effect  the  value 
of  possible  total  profits. 

Next  consider  the  value  of  the  remaining  material  within  the  ultimate 
contour.  This  is  found  to  be  55.279  and  together  with  the  maximum  profit  of 
1*1.427  from  period  1  gives  a  possible  maximum  profit  of  96.706.  This  is 
less  than  the  96. 781  profit  value  achieved  in  the  maximum  total  profit  plan. 
Since  96.706  is  an  upper  bound  on  the  total  maximum  profit  for  the  sequential 
suboptimization  schedule  it  is  clear  that  in  this  case  such  an  approach  would 
be  costly. 

The  last  example  has  illustrated  the  possible  profit  effect  of  sequ¬ 
ential  sub-optimization.  Even  when  discounting  we  cannot  be  certain  of  at¬ 
taining  the  goal  of  maximum  total  profit  unless  consideration  is  given  the 
mining  system  over  the  entire  planning  horizon. 

5.3  MANAGEMENT  CONSIDERATIONS 

The  model  and  techniques  developed  in  this  thesis  provide  management 
with  a  new  decision  tool.  Besides  yielding  the  optimal  production  schedule 
and  an  evaluation  of  the  mineral  deposit,  policies,  technological  changes 
and  system  design  features  may  be  evaluated  as  demonstrated  by  the  numerical 
examples  of  the  previous  section. 
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The  cost  and  effect  of  policies  can  be  readily  obtained  by  inserting 
constraints  describing  such  policies.  Sensitivity  analysis  or  parametric 
analysis  will  aid  the  determination  of  the  effect  of  constraints  on  profit. 
For  example,  assume  for  one  reason  or  another  a  limit  has  been  imposed  on 
capital  expenditures  of  a  certain  type,  say  electric  shovels,  thus  limiting 
available  shovel  hours  per  period.  By  examining  the  multiplier  or  price 
corresponding  to  this  restriction,  the  increase  or  decrease  in  profit  re¬ 
lated  to  a  variation  of  such  a  constraint  could  be  indicated. 

Such  an  analysis  could  be  helpful  in  determining  a  capital  investment 
plan.  If  policy  constraints  are  not  tight,  they  are  not  effectual  in 
governing  the  plan  and  hence  need  not  concern  management  or  operators  until 
the  conditions  change.  Among  the  many  management  considerations  which  may 
be  studied  are  orderly  depletion  policies,  equipment  relocation,  marketing 
and  labor  agreements  which  effect  the  operation. 

Technological  changes  such  as  those  which  affect  the  method  and  results 
of  mining,  concentrating  or  refining  can  be  evaluated  by  considering  them 
as  a  new  activity  or  possibly  only  effecting  the  profit  coefficients  of 
present  activities.  For  example,  if  a  new  method  of  concentrating  is  de¬ 
veloped,  this  could  be  considered  as  adding  a  new  element  to  each  as 

discussed  in  Chapter  1.  Thus  using  the  present  prices,  the  economic  value 
of  the  new  technology  could  be  evaluated. 

Sys  tern  design  features  such  as  mining  methods,  concentrator  design, 
haulage  systems  and  location  of  surface  facilities  may  be  evaluated 
through  use  of  the  proposed  model.  The  design  variations  could  be  incor¬ 
porated  into  the  original  model  by  introducing  new  activities  or  evaluated 
as  separate  systems.  The  location  of  surface  facilities  such  as  waste 
dumps,  concentrator  and  other  shops  can  be  studied  by  employing  the  model. 
Problems  of  relocating  existing  facilities,  which  may  for  some  reason  be 
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within  the  area  of  the  possible  mineable  reserve,  can  be  solved  by  using 
a  dummy  block  with  a  cost  equal  to  the  relocating  cost  [10],  The  best  time 
for  relocating  would  also  be  a  product  of  the  model  with  this  technique. 

Since  the  coefficients  of  the  model  are  based  on  drill  hole  samples, 
their  estimates  may  not  always  be  precisely  accurate.  This  fact  makes 
sensitivity  analysis  [5,  21]  extremely  important  [12].  Many  of  the  co¬ 
efficients  may  be  the  result  of  a  statistical  analysis,  thus  using  sen¬ 
sitivity  analysis  in  conjunction  with  the  confidence  intervals  for  these 
coefficients  could  aid  in  determining  where  expenditures  should  be  made 
to  obtain  better  estimates.  Also,  sensitivity  analysis  is  important  due 
to  the  variability  that  occurs  in  a  mining  operation.  Some  of  this  is  con¬ 
trollable  while  in  other  cases  uncontrollable  c i rcunv, tances  may  cause  dis¬ 
ruptions  in  the  optimal  plan.  Knowing  the  effects  of  variability  through 
sensitivity  analysis  could  aid  greatly  once  a  disruption  occurs  in  the 
plan  of  operation. 

The  great  effect  mine  planning  has  on  corporate  profits  certainly 
justifies  a  scientific  approach  to  the  problem.  The  cost  associated  with 
a  poor  mining  plan  may  far  exceed  the  cost  of  development  and  utilization 
of  new  scientific  decision  techniques  which  greatly  improve  mine  planning 
and  production  scheduling.  Even  though  the  trial  and  error  techniques  which 
have  been  developed  recently  have  improved  mine  planning,  they  still  contain 
many  of  the  inherent  erroneous  concepts  of  the  traditional  methods  as  has 
been  pointed  out  in  this  work.  Often  with  the  trial  and  error  methods, 
many  of  the  important  factors  influencing  profit  are  neglected  and  there 
is  always  that  element  of  doubt  concerning  how  the  chosen  plan  rates  with 
regard  to  the  "best  possible." 


As  a  readily  adaptable  suggestion  to  those  who  are  presently  utilizing 
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extensive  trial  and  error  methods  we  present  the  following.  Formulate  the 
system  constraints  as  a  linear  programming  master  problem  and  use  present 
techniques  of  selecting  a  plan  to  solve  the  subproblems.  The  master  prob¬ 
lem  will  select  and  combine  the  trial  plans  into  a  feasible  schedule  and 
the  pricing  mechanism  will  guide  the  trial  plan  selection  towards  the  op¬ 
timal.  Although  not  as  desirable  as  the  methods  developed  in  this  study, 
such  a  scheme  should  improve  the  trial  and  error  method  with  very  little 
additional  work. 

The  model  and  the  method  of  solution  developed  through  this  research 
have  great  potential  application.  Their  flexibility  greatly  improves  present 
methods  of  mine  planning.  The  profit  derived  from  the  system  approach  to 
the  cut-off  alone  may  pay  for  the  cost  of  practical  implementation.  It 
Is  hoped  that  the  tools  provided  by  these  techniques  will  be  a  benefit  to 
engineers,  operators,  mine  and  corporate  management  through  increased  ef¬ 
ficiency  and  effectiveness  in  planning,  scheduling,  evaluating  and  managing 
their  operations. 
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APPENDIX  A 


NETWORK  DEFINITIONS 

The  following  definitions  of  network  or  graph  terminology  are  an  aid  to 
the  interpretation  of  the  development  in  the  main  text. 

1.  Graph  or  Network  (Directed)  the  system,  E  *  (x,A)  formed  by  a 
set  of  elements  x^  c  x  called  nodes,  and  a  set  of  ordered  pairs 
(x^,Xj)  e  A  called  arcs. 

2.  Sub  graph  of  (x,A)  -  a  graph  (Y,A)  where  Y  c  x  and 

Ay  «  {  (x^Xj)/(XjXj)  e  A,x^,Xj  c  Y}  in  other  words  it  consists 
of  a  subset,  Y  of  x  and  all  the  arcs  in  A  which  connect  tha 
nodes  in  Y  . 

3.  Partial  graph  of  (x,A)  -  a  graph  (x,Ap)  formed  by  the  set  x 

and  the  set  of  arcs  Ap  c  A  ;  i.e.,  a  partial  graph  contains  all 

the  nodes  of  the  original  graph  and  a  partial  set  of  arcs. 

A.  Partial  subgraph  of  (x,A)  -  a  partial  graph  of  a  subgraph  of  the 
form  (Y,A  )  ,  where  Y  c  x  and  A  C  A  c  A  . 

yp  yp  y 

Remark: 

An  arc  (x. ,x4)  or  (l,j)  will  always  be  assumed  to  be  directed  from 

*  J 

or'  i  to  Xj  or  j  . 

5.  Edge  [x^.Xj]  of  the  graph  (x,A)  -  a  pair  of  nodes  [x^,Xj]  is 

an  edge  if  (x^.x^)  c  A  or  (x^.x^  e  A  . 

6.  Adjacent  nodes  -  two  distinct  nodes  joined  by  an  arc. 

7.  Adjacent  arcs  (edges)  -  two  arcs  (edges)  with  a  common  node. 

8.  Path  -  a  sequence  of  arcs  {a  ,  a^  ,  a^  ,  a^)  such  that  the 


terminal  node  of  each  arc  coincides  with  the  initial  node  o?  the 


A.  2 


succeeding  arc;  i.e.,  all  arcs  in  the  path  have  a  similar  orienta¬ 
tion. 


9. 


10. 


11. 


A  path  is  compound  if  it  uses  the  same  arc  more  than  once. 

A  path  is  ele.  jntary  if  it  does  not  pass  through  the  same  node 
more  than  once. 

Circuit  -  a  path  {a,  . aR}  such  that  the  initial  node  of  a^ 
and  the  terminal  node  of  a  coincide. 


n 


Strongly  connected  graph  -  a  graph  in  which  for  any  two  nodes 


and  Xj  ,  x^  i  there  exists  a  path  connecting  x^  with 

XJ  ‘ 

Accessible  sets  ~  when  x^  and  x^  are  nodes  such  that  there  is  a 
path  from  xA  to  x^  ,  then  x^  is  accessible  from  xA  written 
as  x^  >  Xj  .  The  accessible  set  DCx^)  is  the  set  of  all  nodes 


accessible  from  x 


i  ’ 


Note: 

By  definition  xA  {  DCx^)  also  define  »  DCx^  U  x^ 

12.  Chain  -  a  sequence  of  edges  . *R)  *v»ch  that  each 

of  the  intermediate  edges  l ^  is  attached  to  at  one  of  its 

extremities  and  to  ^  at  the  other. 

13.  Cycle  -  a  chain  beginning  and  ending  at  the  same  node. 

14.  Connected  graph  -  a  graph  in  which  for  any  pair  of  nodes  x^  and 
Xj  there  exists  a  chain  connecting  them. 

13.  Antisymmetric  graph  *  a  graph  (x,A)  such  that 
(x^  »Xj )  c  A  o  (Xj,x1)  {  A  . 

16.  Incidence  -  an  arc  whose  Initial  node  is  x  is  said  to  be  "incident 


out  from"  x  -  an  arc  whose  terminal  node  is  x  is  said 


to  be  "incident  into"  x  . 

given  the  set  of  nodes  x  ,  then  the  arc  a^  ■  (x,y)  . 
is  incident  into  x  if  y  e  x  ,  and  x  (  x  . 
and  incident  out  from  x  if  x  c  X  and  y  {  X  , 

m 

The  set  of  arcs  Incident  into  x  is  denoted  as  A(x)  and 

+ 

the  set  incident  out  from  x  as  A(x)  . 

-  + 

17.  Co-boundary  -  the  nonempty  set  of  arcs  A(x)  U  A(x)  . 

18.  Co-circuit  -  a  co-boundary  in  which  all  arcs  are  In  the  same  set, 

+  — 
either  in  A(x)  or  all  in  A(x)  . 

19.  Tree  -  a  connected  graph  which  contains  no  cycleal  It  is  connected 
and  has  n-nodes  and  n  -  1  arcs. 

20.  Node-arc  incidence  matrix  (E)  -  given  a  graph  or  network  {x,A} 
the  node-arc  Incidence  matrix  is  formed  by: 

1.  Listing  all  nodes  vertically  z  row  index 

2.  Listing  all  arc  horlzontlally  =  column  index 

3.  The  elements  of  E  are  formed  by 

.  +  1  if  arc  j  is  incident  out  of  node  i 

■  -  1  if  arc  j  is  incident  into  node  i 
0  otherwise 


For  example  the  node-arc  incidence  matrix  for  the  network  in 
Figure  A.l  is 
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